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Abstract 

Commutative complex numbers of the form u = x + ay + j3z + yt in 4 dimensions are 
studied, the variables x, y, z and t being real numbers. Four distinct types of multiplication 
rules for the complex bases a,f3 and 7 are investigated, which correspond to hypercomplex 
entities called in this paper circular fourcomplex numbers, hyperbolic fourcomplex numbers, 
planar fourcomplex numbers, and polar fourcompex numbers. Exponential and trigonometric 
forms for the fourcomplex numbers are given in all these cases. Expressions are given for the 
elementary functions of the fourcomplex variables mentioned above. Relations of equality 
exist between the partial derivatives of the real components of the functions of fourcomplex 
variables. The integral of a fourcomplex function between two points is independent of the 
path connecting the points. The concepts of poles and residues can be introduced for the 
circular, planar, and polar fourcomplex numbers, for which the exponential forms depend on 
cyclic variables. A hypercomplex polynomial can be written as a product of linear factors 
for circular and planar fourcomplex numbers, and as a product of linear or quadratic factors 
for the hyperbolic and polar fourcomplex numbers. 
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1 Introduction 



A regular, two-dimensional complex number x + iy can be represented geometrically by 
the modulus p = (x 2 + y 2 ) 1 ^ 2 and by the polar angle 9 = arctan(y/x). The modulus p is 
multiplicative and the polar angle 9 is additive upon the multiplication of ordinary complex 
numbers. 

The quaternions of Hamilton are a system of hypercomplex numbers defined in four 
dimensions, the multiplication being a noncommutative operation, fij and many other hy- 
percomplex systems are possible, [||-[|| but these hypercomplex systems do not have all the 
required properties of regular, two-dimensional complex numbers which rendered possible 
the development of the theory of functions of a complex variable. 

This paper belongs to a series of studies on commutative complex numbers in n dimen- 
sions. || Systems of hypercomplex numbers in 4 dimensions of the form u = x + ay + (3z + jt 
are described in this work, where the variables x, y, z and t are real numbers, for which 
the multiplication is both associative and commutative. The product of two fourcomplex 
numbers is equal to zero if both numbers are equal to zero, or if the numbers belong to 
certain four-dimensional hyperplanes as discussed further in this work. The fourcomplex 
numbers are rich enough in properties such that an exponential and trigonometric forms 
exist and the concepts of analytic fourcomplex function, contour integration and residue 
can be defined. Expressions are given for the elementary functions of fourcomplex vari- 
able. The functions f(u) of fourcomplex variable defined by power series, have derivatives 
lim u ^ Uo [f(u) — f(uo)]/(u—uo) independent of the direction of approach of u to uq. If the four- 
complex function f(u) of the fourcomplex variable u is written in terms of the real functions 
P(x, y, z, t), Q(x, y, z, t),R(x, y, z, t), S(x, y, z, t), then relations of equality exist between par- 
tial derivatives of the functions P, Q, R, S. The integral fj± f(u)du of a fourcomplex function 
between two points A, B is independent of the path connecting A, B. 

Four distinct types of hypercomplex numbers are studied, as discussed further. In Sec. 
II, the multiplication rules for the complex units a,/3 and 7 are a 2 = —1, [5 2 = —1, 7 2 = 
l,a/3 = Pa = —7, 07 = ja = /3, (3~f = 7/3 = a. The exponential form of a fourcomplex 
number is x + ay + (5z + 7* = pexp [7 In tan?/; + (l/2)a(4> + \) + (l/2)/3(0 — x)] > where 
the amplitude is p 4 = [(x + t) 2 + (y + z) 2 ] [(x — t) 2 + (y — z) 2 ] , (j), x are azimuthal angles, 
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< <fi < 2tt, < x < 2-7T, and ^ is a planar angle, < ip < tt/2. The trigonometric 
form of a fourcomplex number is x + ay + /3z + jt = d[cos(ip — 7r/4) + 7sin(?/> — vr/4)] 
exp [(l/2)a(0 + x) + (l/2)/?(0 - x)], where d 2 = x 2 + y 2 + z 2 + t 2 . The amplitude p and 
tan ^ are multiplicative and the angles 4>, \ are additive upon the multiplication of fourcom- 
plex numbers. Since there are two cyclic variables, <p and X-> these fourcomplex numbers 
are called circular fourcomplex numbers. If f(u) is an analytic fourcomplex function, then 
§ v f{u)du/(u - n ) = vr[(a + 0) int(n ^,r^) + (a - (3) int(u 0r £ , T tC )} f(u ), where the 
functional int takes the values or 1 depending on the relation between the projections of 
the point uq and of the curve T on certain planes. A polynomial can be written as a product 
of linear or quadratic factors, although the factorization may not be unique. 

In Sec. Ill, the multiplication rules for the complex units a, (3 and 7 are a 2 = 1, (3 2 = 
1, 7 2 = 1, a(3 = (3a = 7, 07 = 7a = j3, (3^ = ^(3 = a. The exponential form of a fourcomplex 
number, which can be defined for s = x+y+z+t > 0, s' = x — y+z — t > 0, s" = x+y — z — t > 
0, s'" = x — y — z + t > 0, is x + ay + f3z + = p exp(cn/i + (3z\ +7ti), where the amplitude is 
p = (ss 1 V s'") 1 / 4 and the arguments are y\ = (1/4) h\(ss" / s' s'"), z\ = (1/4) ln(ss"' /s's"), t\ = 
(1/4) ln(ss' /s"s"'). Since there is no cyclic variable, these fourcomplex numbers are called 
hyperbolic fourcomplex numbers. The amplitude p is multiplicative and the arguments 
yi, z\,t\ are additive upon the multiplication of fourcomplex numbers. A polynomial can be 
written as a product of linear or quadratic factors, although the factorization may not be 
unique. 

In Sec. IV, the multiplication rules for the complex units a, [3 and 7 are a 2 = (3, [3 2 = 
— 1, 7 2 = —(3,a(3 = (3a = 7, cry = 7a = —1, /?7 = 7/? = —a . The exponential 
function of a fourcomplex number can be expanded in terms of the four-dimensional co- 
sexponential functions fw{x) = 1 — x 4 /4! + x 8 /8\ — • • • , f&i(x) = x — x 5 /5! + x 9 /9! — 
■ • ■ , / 42 (x) = x 2 /2\ - x 6 /6! + x 10 /10! , / 43 (x) = x 3 /3! - x 7 /7! + x n /W. . Ex- 
pressions are obtained for the four-dimensional cosexponential functions in terms of ele- 
mentary functions. The exponential form of a fourcomplex number is x + ay + (3z + jt = 
p exp {(1/2) (a - 7) In tan ^ + (1/2) [0 + (a + 7)/^]^ - (1/2) [(3 - (a + 7)/V^]x}. where the 
amplitude is p 4 = j [x + (y - t)/y/2] 2 + [z + (y + t)/y/5\ 2 J j [x - (y - t)/y/2~\ 2 + [z - (y + t)/v^] 2 J, 
0, x are azimuthal angles, < e6 < 2-7T, < x < 2-7T, and ip is a planar angle, < V < 7r /2- 
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The trigonometric form of a fourcomplex number is x + ay + [3z + = d [cos (ip — vr /4) 
+(1/V2)(a - 7) sin - tt/4)] exp {(l/2)[/3 + (a + 7)/^ - (1/2) [/3 - (a + 7)/V%}, where 
d 2 = x 2 + y 2 + z 2 + t 2 . The amplitude p and tan^ are multiplicative and the angles 4>, \ 
are additive upon the multiplication of fourcomplex numbers. There are two cyclic vari- 
ables, 4> and x> so that these fourcomplex numbers are also of a circular type. In order 
to distinguish them from the circular hypercomplex numbers, these are called planar four- 
complex numbers. If f(u) is an analytic fourcomplex function, then § r f(u)du/(u — uq) = 
7r \B+ (a + 7)/\/2j int(no^,r^) + ( f3 - (a + rf/VV int(n 0T ^, T Tf ) f(u ) , where the 
functional int takes the values or 1 depending on the relation between the projections 
of the point no and of the curve T on certain planes. A polynomial can be written as a 
product of linear or quadratic factors, although the factorization may not be unique. The 
fourcomplex numbers described in this section are a particular case for n = 4 of the planar 
hypercomplex numbers in n dimensions. ||,|| 

In Sec. V, the multiplication rules for the complex units a, (3 and 7 are a 2 = (3, (3 2 = 
1, 7 2 = (3,a(3 = Pa = 7, 07 = 7a = 1, /?7 = 7/? = a. The product of two fourcomplex 
numbers is equal to zero if both numbers are equal to zero, or if the numbers belong to cer- 
tain four-dimensional hyperplanes described further in this work. The exponential function 
of a fourcomplex number can be expanded in terms of the four-dimensional cosexponential 

functions g4 (x) = 1 + x 4 /4! + x 8 /8! H , 541 (x) = x + x 5 /5l + x 9 /9! H , g 42 (x) = x 2 /2! + 

x 6 /6! + x 10 /10! + • • • , <?43(x) = x 3 /3! + x 7 /7! + x /ll! + • • •. Addition theorems and other 
relations are obtained for these four-dimensional cosexponential functions. The exponential 
form of a fourcomplex number, which can be defined for x+y+z+t > 0, x—y+z—t > 0, is u = 



pexp 



where 



(1/4) (a + p + 7) ln^/tan 6+) - (1/4) (a - P + 7) ln(>/2/ tan 0_) + (a - 7 )0/2 
p = {p+p-) 1/2 , /4 = (x-z) 2 +{y-t) 2 , /1 2 = {x+z) 2 -{y+t) 2 , e + = (l+a+/3+ 7 )/4, e_ = (1- 
a+P~ 7)/4, ei = (1 — P)/2,e\ = (a — 7)/2, the polar angles are tan# + = \f2p + /v +: t&\i6- = 
\/2/i+/i>-, < 9+ < ir, < 6- < 7r, and the azimuthal angle is defined by the relations 
x — y = p+ cos 0, z — t = p + sin (j), < <p < 2ir. The trigonometric form of the fourcomplex 
number u is u = (1 + 1/tan 2 6 + + 1/tan 2 6L)~ 1/2 |ei + e + \/2/tan0+ + e_\/2/tan6L j 
exp [ei</)]. The amplitude p and tan 6+/\/2, tan #_/\/2, are multiplicative, and the azimuthal 
angle </> is additive upon the multiplication of fourcomplex numbers. There is only one 



4 



cyclic variable, (f), and there are two axes v+,V— which play an important role in the 
description of these numbers, so that these hypercomplex numbers are called polar four- 
complex numbers. If f(u) is an analytic fourcomplex function, then § r f(u)du/(u — uq) = 
7r(/3 — 7) mt(uQ^ v , T^ v )/(uq) , where the functional int takes the values or 1 depending on 
the relation between the projections of the point uq and of the curve T on certain planes. A 
polynomial can be written as a product of linear or quadratic factors, although the factoriza- 
tion may not be unique. The fourcomplex numbers described in this section are a particular 
case for n = 4 of the polar hypercomplex numbers in n dimensions. 

2 Circular Complex Numbers in Four Dimensions 
2.1 Operations with circular fourcomplex numbers 

A circular fourcomplex number is determined by its four components (x,y,z,t). The sum 
of the circular fourcomplex numbers (x,y,z,t) and (x' , y' , z', t') is the circular fourcomplex 
number (x + x',y + y',z + z',t + t'). The product of the circular fourcomplex numbers 
(x,y,z,t) and (x' , y', z', t') is defined in this work to be the circular fourcomplex number 
(xx' — yy' — zz' + tt' , xy' + yx' + zt' + tz' , xz' + zx' + yt' + ty\ xt' + tx' — yz' — zy'). 

Circular fourcomplex numbers and their operations can be represented by writing the 
circular fourcomplex number (x, y, z,t) asu = x + ay + (3z + jt, where a, (3 and 7 are bases 
for which the multiplication rules are 

a 2 = —1, 1 = —1, 7 2 = 1, a/3 = j3a = —7, 07 = 7a = (3, /?7 = 7/3 = a. (1) 

Two circular fourcomplex numbers u = x + ay + (3z + "ft, v! = x' + ay' + f3z' + yt' are equal, 
u = u', if and only if x = x' , y = y' , z = z' , t = t' . If u = x+ay+(3z+yt, u' = x' '+ay' '+{3z' '+ r yt' 
are circular fourcomplex numbers, the sum u + u' and the product uv! defined above can be 
obtained by applying the usual algebraic rules to the sum (x + ay + (3z + yt) + (x' + ay' + 
j3 z' + 7i') and to the product (x + ay + (3z + yt){x' + ay' + f3z' + yt 1 ), and grouping of the 
resulting terms, 

u + u ' = x + x' + a(y + y') + (3{z + z) + y(t + t'), (2) 
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uu = xx — yy — zz + tt + a{xy + yx' + zt' + tz ) + /3(xz' + zx + yi + ty') 

+7(xt' + ix' — yz' — zy') (3) 

If u,u',u" are circular fourcomplex numbers, the multiplication is associative 

(uu')u" = u(u'u") (4) 
and commutative 

uu' = u'u, (5) 

as can be checked through direct calculation. The circular fourcomplex zero is + a ■ + 
j3 • + 7 • 0, denoted simply 0, and the circular fourcomplex unity is l + a- + /3-0 + 7-0, 
denoted simply 1. 

The inverse of the circular fourcomplex number u = x + ay + (3z + jt is a circular 
fourcomplex number u' = x' + ay' + + jf having the property that 

uu' = 1. (6) 

Written on components, the condition, Eq. (^), is 

xx' — yy' — zz' + tt' = 1, 
yx' + xy' + tz' + zt' = 0, 
zx' + ty' + xz' + yt! = 0, 
ix' — zy' — yz' + xi' = 0. 
The system (^) has the solution 

, _ x(x 2 + y 2 + z 2 - t 2 ) - 2yzt 

, _ y(~x 2 -y 2 + z 2 - t 2 ) + 2xzt 

y a i v y j 

p 

, _ z(-x 2 + y 2 -z 2 -t 2 ) + 2xyt 



(7) 



, _ i(-x 2 + y 2 + z J + r) - 2xyz 

t= 7 : 



ii 



provided that p ^ 0, where 

4 4, 4, 4 , .4 , 2 2 , 22 2,2 2 2, 2,2 , 2,2n o j. /in\ 

p = x +y + z +t + 2(x y + x z — x i — y z + y t + z t ) — Sxyzt. (12) 



The quantity p will be called amplitude of the circular fourcomplex number x + ay + (3z + yt. 
Since 

P 4 = P 2 + P 2 -, (13) 

where 

p\ = (x + tf + (y + zf, pi = {x- t) 2 + {y- z)\ (14) 
a circular fourcomplex number u = x + ay + (3z + jt has an inverse, unless 

x + t = 0, y + z = 0, (15) 

or 

x-t = 0, y-z = 0. (16) 



Because of conditions (|l5|)-(16) these 2-dimensional surfaces will be called nodal hyper- 
planes. It can be shown that if uvf = then either u = 0, or a' = 0, or one of the 
circular fourcomplex numbers is of the form x + ay + /3y + 72; and the other of the form 
x' + ay' — 0y' — 7a/. 

2.2 Geometric representation of circular fourcomplex num- 
bers 

The circular fourcomplex number x + ay + (3z + jt can be represented by the point A of 
coordinates (x, y, z, t). If O is the origin of the four-dimensional space x, y, z, t, the distance 
from A to the origin O can be taken as 

d 2 = x 2 + y 2 + z 2 + t 2 . (17) 

The distance d will be called modulus of the circular fourcomplex number x + ay + f3z + jt, 
d = \u\. The orientation in the four-dimensional space of the line OA can be specified with 
the aid of three angles </>, x> Y> defined with respect to the rotated system of axes 
x + i x-t y + z y - z 

t = l^> T = l^> v = ^T> c = ^T- (18) 

The variables £, t>, r, £ will be called canonical circular fourcomplex variables. The use of the 
rotated axes £, u, r, £ for the definition of the angles </>, "0 is convenient for the expression 



of the circular fourcomplex numbers in exponential and trigonometric forms, as it will be 
discussed further. The angle <p is the angle between the projection of A in the plane £, v 
and the 0£ axis, < (ft < 2ir, x is the angle between the projection of A in the plane r, ( 
and the Or axis, < x < 2vr, and ip is the angle between the line OA and the plane tO(, 
< ip < vr/2, as shown in Fig. 1. The angles (ft an d X will be called azimuthal angles, the 
angle ip will be called planar angle . The fact that < ip < it/ 2 means that ip has the same 
sign on both faces of the two-dimensional hyperplane vOQ. The components of the point A 
in terms of the distance d and the angles (ft, %, ip are thus 
x + t 



V2 
x — t 

y + z 

y-z 



= dcoscj)sm-ip, (19) 

= d cos x cos ip, (20) 

= <i sin sin ip, (21) 

= dsinxcosip. (22) 



V2 

It can be checked that p + = \/2dsmtp, p- = \/2dcosip. The coordinates x,y,z,t in terms of 
the variables d, (ft, x, ip are 

x = -^= (cos ip cos x + sin ip cos (ft) , (23) 

y = -^= (cos ip sin % + sin ip sin 0) , (24) 

z = —t=(— cos V'sinx + sin ?/> sin 0) , (25) 
V2 

t = — /=(— cos ^ cos x + sin?/> cos </>). (26) 
V2 

The angles (ft, x, ip can be expressed in terms of the coordinates x, y, z, t as 

sin(ft= (y + z)/p + , cos<p= (x + t)/p + , (27) 

sin x = {y~ z)/p- , cos x = (x ~ t)/p- , (28) 
tanip = p + /p_. (29) 
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The nodal hyperplanes are £Ov , for which t = 0, £ = 0, and tOQ, for which £ = 0, v = 0. 
For points in the nodal hyperplane £Ov the planar angle is ift = vr/2, for points in the nodal 
hyperplane tOQ the planar angle is ift = 0. 

It can be shown that if % = x\ + ayi + /Jzi + 7*1, v>2 = X2 + «2/2 + j3z 2 + 7*2 are circular 
fourcomplex numbers of amplitudes and angles pi, </>i, xi, i/>i and respectively P2,4>2,X2,'>p2, 
then the amplitude p and the angles (ft, x, ip of the product circular fourcomplex number u\U2 
are 



The relations are consequences of the definitions ©-© and of the 



P = P1P2, 



(30) 



(ft = <fti + (ft: 



'2, X = Xi + X2, tanip = tan -01 tan ^2- 



(31) 



identities 



[(xix 2 - 2/12/2 - ziz 2 + ht 2 ) + (xit 2 + hx 2 - yiz 2 - ziy 2 )f 

+ [{x\y 2 + y\X2 + Z X t2 + hz 2 ) + (X\Z2 + Z\X 2 + 2/1*2 + *l2/2)] 2 

= \{x x + h) 2 + (yi + zi) 2 ][(x 2 + t 2 ) 2 + (2/2 + z 2 )\ 



(32) 




(33) 



(x x x 2 - J/13/2 - Zi*2 + M2) + (a?i*2 + *i^2 - 2/1*2 - ^12/2) 



(xi + ii)(ar 2 + h) - (2/1 + zi)(y 2 + 22) 



(34) 



{x x x 2 - 2/12/2 - *i*2 + ht 2 ) - (xit 2 + t\x 2 - y\z 2 - ziy 2 ) 



{x\ - ti)(x 2 - t 2 ) - (2/1 - zi)(y 2 - z 2 ) 



(35) 



(xiV2 + V1X2 + zit 2 + tiz 2 ) + (xiz 2 + zix 2 + 2/1*2 + £12/2) 



(2/1 + 2l)(iC2 + + (2/2 + Z 2 ){X2 + t 2 ) 



(36) 



(#12/2 + 2/1X2 + 21*2 + *1*2) - (^1*2 + ZiX 2 + 2/1*2 + *l2/2) 

= (2/1 - zi)(x 2 - h) + (2/2 - 2 2 )(x 2 - t 2 ). 
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(37) 



The identities (p2j) and (p3f) can also be written as 



p + = P1+P2+, 



P- = P1-P2- 



(38) 
(39) 



where 



+ tj) 2 + (y j+Zj ) 2 , p)_ 



(40) 



The fact that the amplitude of the product is equal to the product of the amplitudes, as 



written in Eq. (30), can be demonstrated also by using a representation of the multiplication 
of the circular fourcomplex numbers by matrices, in which the circular fourcomplex number 
u = x + ay + j3z + is represented by the matrix 



A 



( + \ 

x y z t 

—y x t —z 

—z t x —y 



(41) 



\ t z y x J 

The product u = x + ay + fiz + jt of the circular fourcomplex numbers u\ = x\ + ay\ + 
/3z% + 7*1, U2 = X2 + ay2 + (3z2 + 7*21 can be represented by the matrix multiplication 



A = A X A 2 . 

It can be checked that the determinant det(j4) of the matrix A is 
deU = p 4 . 



(42) 



(43) 



The identity ([K]) is then a consequence of the fact the determinant of the product of matrices 
is equal to the product of the determinants of the factor matrices. 



2.3 The exponential and trigonometric forms of circular four- 
complex numbers 

The exponential function of the fourcomplex variable u is defined by the series 



exp u = 1 + u + u 2 /2l + n 3 /3! + 



(44) 
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It can be checked by direct multiplication of the series that 

exp(u + u ) = exp u ■ exp u . (45) 

These relations have the same form for all hypercomplex systems discussed in this work. 

If u = x + ay + (3z + jt, then exp u can be calculated as exp u = exp x • exp (ay) • exjp(ftz) ■ 
exp(7i). According to Eq. ([[]), 

a =(-1) ,a x =(-1) a,p =(-1) ,p =(-1) P,7 = 1> ( 46 ) 
where m is a natural number, so that exp(ay), exp(/?z) and exp(7t) can be written as 

exp(ay) = cos y + a sin y, exp(/?z) = cos z + /? sin z, (47) 

and 

exp(7£) = cosh t + 7 sinh i. (48) 

From Eqs. (|47|)-([48|) it can be inferred that 

(cos y + a sin y) m = cos my + a sin my, 

(cos 2; + /3sinz) m = cos mz + /3 sin mz, (49) 
(cosh i + 7 sinh t) m = cosh mi + 7 sinh mi. 

Any circular fourcomplex number u = x + ay + @z + jt can be writen in the form 

x + ay + /?z + jt = e Xl+ayi +/3zi + 7 * 1 . (50) 

The expressions of x±, y±, z\,t\ as functions of x, y, z, t can be obtained by developing e ayi , e^ Zl 
and e 7 * 1 with the aid of Eqs. (|47|) and (fig), by multiplying these expressions and separating 
the hypercomplex components, 

x = e Xl (cosyi coszi coshti — sinyi sinzi sinhii), (51) 

y = e xi (sin y% cos z\ cosh i 1 + cos yi sin zi sinh t\) , (52) 

z = e xi (cos yi sin z\ cosh ti + sin yi cos z\ sinh i 1 ) , (53) 

t = e Xl (— sinyi sinzi coshti + cosyi coszi sinhii), (54) 
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From Eqs. (|5l|)-(p4|) it can be shown by direct calculation that 

x 2 + y 2 + z 2 + t 2 = e 2xi cosh2ti, (55) 

2(xt + yz) = e 2xi sinh 2tx, (56) 

so that 

By comparing the expression in the right-hand side of Eq. ( |57D with the expression of p, Eq. 
fll3|), it can be seen that 

e Xl = p. (58) 

The variable t\ is then given by 

d 2 , 2(x£ + yz) . , 

cosh2ti = sinh2ti = 5 . (59) 

p z p z 

From the fact that p 4 = d 4 — A(xt + yz) 2 it follows that d 2 /p 2 > 1, so that Eq. fl59|) has 

always a real solution, and t\ = for xt + yz = 0. It can be shown similarly that 

x 2 - y 2 + z 2 - t 2 . 2(xy - g) 

cos2yi = 5 , sin2yi = 3 , (60) 

p z p z 

x 2 + y 2_ z 2_ t 2 2(xz - yt) 

coszzi = 7i , smzzi = 7. . (61) 

p z p A 

It can be shown that (x 2 — y 2 + z 2 — t 2 ) 2 < p 4 , the equality taking place for xy = zt, and 
(x 2 + y 2 — z 2 — t 2 ) 2 < p 4 , the equality taking place for xz = yt, so that Eqs. ( |60|) and Eqs. 
(|6l|) have always real solutions. 
The variables 

1 2(xy - zt) 1 2(xz -yt) 1 , 2(xt + yz) 

yi = - arcsm = , z\ = - arcsm 5 , t\ = -argsmh 7. (62) 

2 p z 2 p z 2 p z 

are additive upon the multiplication of circular fourcomplex numbers, as can be seen from 
the identities 

(xx' - yy' - zz' + tt')(xy' + yx' + zt' + tz') 

— (xz' + zx + yt' + ty')(xt' — yz' — zy + tx ) 

= (xy - zt)(x' 2 - y' 2 + z' 2 - t' 2 ) + (x 2 - y 2 + z 2 - t 2 )(x'y' - z't'), (63) 
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{xx — yy' — zz + tt ){xz + zx + yt + ty ) 

— {xy + yx' + zt' + tz')(xt' — yz — zy + tx) 

I ±\( 12 , 12 12 j_i2\ , / 2 i 2 2 j.2w / / 

= (xz - yt)(x +y - z -t ) + (x +y - z - t )(x z - y t ), 



(64) 



(xx' — yy' — zz + tt )(xt — yz — zy + tx) 

+ (xy' + yx' + zt' + tz')(xz + zx + yt' + ty') 

= (xt + yz)(x' 2 + y' 2 + z' 2 + t' 2 ) + (x 2 + y 2 + z 2 + i 2 )(xY + y'z'). (65) 

The expressions appearing in Eqs. (|5S|)-(l6l"l) can be calculated in terms of the angles 
i, x> ^ with the aid of Eqs. (p3|)-(p6|) as 



1 2(xt + yz) 



d 2 



p 2 sin 2ip ' p 2 



1 



tan 2ip ' 



x 2 — y 2 + z 2 — t 2 

7 



2(xy — zt) 

COS{(f> + X), o = Sm (^ + X), 

p 2 



a; 2 + y 2 - z 2 - t 2 2(xz - yt) . 

2 = cos (^ ~ X), o = sm W ~ X)- 

p z p z 

Then from Eqs. (H)-© and (H)-@ it results that 
Vi = 2 , Zl = 2 , h = -lntam/>, 



(66) 
(67) 
(68) 



so that the circular fourcomplex number u, Eq. (|50|), can be written as 

' 4> -\- x 4> ~ x i 

ii = pexp a — - h/3 — - h 7- hi tan ^ 



(69) 



(70) 



In Eq. (70) the circular fourcomplex number u = x + ay + j3z + jt is written as the product 
of the amplitude p and of an exponential function, and therefore this form of u will be called 
the exponential form of the circular fourcomplex number. It can be checked that 



exp 



exp 



a + [3 



a- 13 



-X 



I-7 1 + 7 a + (3 . 
— r 1 — cos <t> H 8111 ( 



1 + 7 1-7 a- 13 . 
H — cos x H ~ — sin < 



(71) 
(72) 



2 7 V 2 2 2 

which shows that e <ya+ ^ c ^/ 2 and e( Q_ ^ X//2 are periodic functions of <fi and respectively x> with 
period 2-7T. 
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The relations between the variables yi,z±,tx and the angles (j),x,4 ) can De obtained al- 
ternatively by substituting in Eqs. (|5T|)-(|54|) the expression e»i = d/(cosh2t 1 ) 1 /2 ) Eq. <M) 
and summing and subtracting of the relations, 
x + t 



V2 

x — t 
^J2 



dcos(yi + z\) sin(r/ + vr/4), (73) 
dcos(yi — z\) cos(r? + vr/4), (74) 



y ~\~ z 

—j=- = d sin(yi + z\ ) sin(?7 + 7r/4) , (75) 

y z = dsin(yi - z{) cos(r? + 7r/4), (76) 
v2 

where the variable rj is defined by the relations 

coshti sinhti . 

cos 7], , = sin 7/, (77) 



(cosh2ti) 1 /2 " (cosh2ti) 1 /2 

and when — oo < t\ < oo, the range of the variable r] is — 7r/4 < r/ < vr/4. The comparison of 
Eqs. ©-(IID and ©"© shows that 



<t> = Vl + zi, X = V\ - zi, i> = ■n + n/A. (78) 
It can be shown with the aid of Eq. (^8|) that 

ex P Q7 In tan i^j = ( sin 2^)1/2 l cos (^ - T""/ 4 ) + 7 sin(^ - vr/ 4)] . (79) 
The circular fourcomplex number it, Eq. (fFo|), can then be written equivalently as 



it = d{cos(ij) — vr/4) + 7sin('0 — vr/4)} exp 



+ X i-x 
a — h /x 



(80) 



2 ' 2 

In Eq. (^), the circular fourcomplex number u = x + ay + /3z + jt is written as the product 
of the modulus d and of factors depending on the geometric angles eft, x and tp, and this form 
will be called the trigonometric form of the circular fourcomplex number. 

If ui,U2 are circular fourcomplex numbers of moduli and angles d±, 0i, xij "01 and respec- 
tively d2, 4>2: X2, i>2, the product of the planar factors can be calculated to be 

[cos('0i — 7r/4) + 7sin( , 0i — vr/4)] [cos(i/ ; 2 — 7r /4) + 7Sin(^2 — vr/4)] 

= [cos(V>i - i>2) -^cos^i + i/> 2 )]. (81) 
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The right-hand side of Eq. ( |Sl"| ) can be written as 

cos(V>i - 1P2) - 7cos(V'i + 1P2) 

= [2(cos 2 ipi cos 2 V>2 + sin 2 V>i sin 2 V>2)] 1/2 [cos(^ - vr/4) + 7sin(V> - vr/4)], (82) 

where the angle ip, determined by the condition that 

tan(i/> — 7r/4) = — cos(V>i + ^2)/ cos(^i — 1P2) (83) 



is given by tan^ = tan^i tan ip2 1 which is consistent with Eq. (|3l[). It can be checked that 
the modulus d of the product u\U2 is 

d = \f2d\d2 (cos 2 tpi cos 2 ip2 + sin 2 V 7 ! sin 2 ^2) ^ • (84) 



2.4 Elementary functions of a circular fourcomplex variable 

The logarithm u\ of the circular fourcomplex number u, U\ = Inn, can be defined as the 
solution of the equation 

u = e u \ (85) 

written explicitly previously in Eq. (|50|), for U\ as a function of u. Prom Eq. ([7(]) it results 
that 

ln-u = lnp + -7 In tan?/; + a^—^- + [5^—^-. (86) 



It can be inferred from Eqs. (|30| ) an d ( pl|) that 

m(mi') = In u + In u', (87) 

up to multiples of 7r(a + (3) and 7r(a — j3). 

The power function ii n can be defined for real values of n as 

u m = e mlnu . (88) 

The power function is multivalued unless n is an integer. For integer n, it can be inferred 



from Eq. (87) that 



(uuT = u n u' m . 
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If, for example, m = 2, it can be checked with the aid of Eq. ( p0[ ) that Eq. (pq) gives indeed 
(x + ay + (3z + jt) 2 = x 2 - y 2 - z 2 + t 2 + 2a(xy + zi) + 2f3(xz + yt) + 27(2xt - yz). 
The trigonometric functions of the fourcomplex variable u are defined by the series 

cosu = l-u 2 /2\ + u 4 /4\ + --- , (90) 

sin u = u - u 3 /3\ + u 5 /5! + • • • . (91) 

It can be checked by series multiplication that the usual addition theorems hold also for the 
circular fourcomplex numbers u,u', 

cos(u + u ) = cos u cos u — sin u sin u , (92) 

sin(w + n') = sin it cos u + cos u sin u . (93) 

These relations have the same form for all systems of hypercomplex numbers discussed in 
this work. The cosine and sine functions of the hypercomplex variables ay, f3z and jt can be 
expressed as 

cos ay = cosh y, sin ay = a sinh y , (94) 
cos j3y = coshy, sin/Jy = /Jsinhy, (95) 
cos7y = cosy, sin7y = 7siny. (96) 
The cosine and sine functions of a circular fourcomplex number x + ay + (3z + jt can then be 



expressed in terms of elementary functions with the aid of the addition theorems Eqs. (92), 
( |93| ) and of the expressions in Eqs. (|94"|)- (p6|) , 

The hyperbolic functions of the fourcomplex variable u are defined by the series 

cosh u = 1 + u 2 /2\ + n 4 /4! + • • • , (97) 

sinh u = u + ii 3 /3! + n 5 /5! H . (98) 

It can be checked by series multiplication that the usual addition theorems hold also for the 
circular fourcomplex numbers u,u', 

cosh(u + u) = cosh u cosh u + sinh u sinh u , (99) 
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sinh(ii + u ) = sinh u cosh u + cosh u sinh u'. (100) 

These relations have the same form for all systems of hypercomplex numbers discussed in 
this work. The hyperbolic cosine and sine functions of the hypercomplex variables ay, j3z 
and 'ft can be expressed as 

cosh ay = cosy, sinh ay = asiny, (101) 

cosh/3y = cosy, sinh/3y = /3siny, (102) 

cosily = coshy, sinh7y = 7 sinh y. (103) 

The hyperbolic cosine and sine functions of a circular fourcomplex number x + ay + fiz + jt 
can then be expressed in terms of elementary functions with the aid of the addition theorems 
Eqs. (H), ( |T00|) and of the expressions in Eqs. (|IoTl)-( |lO^) . 

2.5 Power series of circular fourcomplex variables 

A circular fourcomplex series is an infinite sum of the form 

a + ai + a 2 H h a n H , (104) 

where the coefficients a n are circular fourcomplex numbers. The convergence of the series 



( 104 ) can be defined in terms of the convergence of its 4 real components. The conver- 
gence of a circular fourcomplex series can however be studied using circular fourcomplex 
variables. The main criterion for absolute convergence remains the comparison theorem, but 
this requires a number of inequalities which will be discussed further. 

The modulus of a circular fourcomplex number u = x + ay + (3z + jt can be defined as 

| u | = ( x 2 +y 2 + z 2 +t 2 ) l/2 ; (1Q5) 

so that, according to Eq. (]T7|), d = \u\. Since \x\ < \u\, |y| < \u\, \z\ < \u\, \t\ < \u\, a property 
of absolute convergence established via a comparison theorem based on the modulus of the 



series (104) will ensure the absolute convergence of each real component of that series. 

The modulus of the sum u\ + u 2 of the circular fourcomplex numbers u\ , u 2 fulfils the 
inequality 

H^il — l u 2|| < |^i + u 2 \ < + \u 2 \. (106) 
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For the product the relation is 

\u\u 2 \ < v2\ui\\u2\, 



(107) 



which replaces the relation of equality extant for regular complex numbers. The equality in 
Eq. ( |1Q7| ) takes place for x\ = t\,y\ = Z\,x<i = *2 , J/2 = %i or x\ = —t\,yi = — Z\,X2 = 
—t2,U2 = —z-2- In Eq. (Bj|), this corresponds to ipi = 0, "02 =0 or -01 = ^/2,ip 2 = 
The modulus of a product, which has the property that < |itiit2|j becomes equal to zero 
for xi = ti, 2/1 = zi,x 2 = -t 2 ,V2 = ~Z2 or x\ = -ii,yi = -zi,«2 = £2,2/2 = ^2, as 



discussed after Eq. (16). In Eq. (84), the latter situation corresponds to ifti =0,-02 = tt/2 
or ipi = 0, -02 = tt/2. 
It can be shown that 



x: 



+ y 2 + z 2 + t 2 <\u 2 \< V2(x 2 + y 2 + z 2 + t 2 



(108) 



The left relation in Eq. (108) becomes an equality, x 2 + y 2 + z 2 + t 2 = \u 2 \, for xt + yz = 0. 
This condition corresponds to tpi = tp 2 = 7r/4 in Eq. fl84]). The inequality in Eq. ( |107| ) 
implies that 



| u '| < 2 ('- 1 )/ 2 |u|'. 
From Eqs. (|l07|) and (BM) it results that 



an 



Ii < 2'/2| r ll..lJ 



am 



A power series of the circular fourcomplex variable u is a series of the form 

.1 i 



ao + aiti + a2n 



+ aiu 



Since 



1=0 



<E 2 ' /2 NM', 

1=0 



a sufficient condition for the absolute convergence of this series is that 
lim ^+ lN < 1. 

l^oo \ai\ 

Thus the series is absolutely convergent for 



\u\ < c, 



(109) 



(110) 



:ni) 



(112) 



(113) 



(114) 
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where 



c = lim p . (115) 

i^oo V2|a/ + i| 



The convergence of the series ( |lll|) can be also studied with the aid of the transformation 
x + ay + (3z + -ft = \/2(ei£ + e x v + e 2 r + e 2 Q, (116) 



where £, v, r, £ have been defined in Eq. (|18|), and 



ei = — g— , ei = — ^— , e 2 = e 2 = . (117) 

The ensemble ei,ei,e 2 ,e 2 will be called the canonical circular fourcomplex base, and Eq. 



(|116|) gives the canonical form of the circular fourcomplex number. It can be checked that 



e? = ei, e\ = -ei, ei&L = ei, e 2 = e 2 , e\ = -e 2 , e 2 e 2 = e 2 , 



eie 2 = 0, eie 2 = 0, eie 2 = 0, e 2 e\ = 0. (118) 



The moduli of the bases in Eq. (|117| ) are 

|eil = 7i' |gll = 7f' |e2l = 7i' |e ~ 2l = 71' (119) 

and it can be checked that 

\x + ay + /3z + jt\ 2 = f + v 2 + T 2 + ( 2 . (120) 

If u = u'u" ', the components £, t>, r, £ are related, according to Eqs. (|34])-(|37|) by 

£ = V2(£e-v'v"),v = V2£v"+v'Z"), r = V2(t't"-C'("), C = V2(t'("+('t"), (121) 

which show that, upon multiplication, the components £, v and r, (" obey, up to a normaliza- 
tion constant, the same rules as the real and imaginary components of usual, two-dimensional 
complex numbers. 



If the coefficients in Eq. (Ill) are 

ai = a ro + ana + Pa i2 + ja t3 , (122) 

and 

An = aio + ai 3 , An = an + a i2 , A i2 = ai - a/ 3) A l2 = an - a t2 , (123) 
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the series (111) can be written as 

oo 

2 l/2 \(eiAn + eii n )(ei£ + e x v) 1 + (e 2 A l2 + e 2 A l2 ){e 2 T + e 2 Q l 



Thus, the series in Eqs. (|111|) and (|124|) are absolutely convergent for 



P+ < ci, p_ < c 2 , 



where 



ci = lim ■ 



1/2 



1/2' C 2 = , lim ' 



^2 + ^2 



1/2 



1/2- 



(124) 



(125) 



(126) 



It can be shown that c = (l/v / 2)min(ci, c 2 ), where min designates the smallest of the 
numbers c\,c 2 . Using the expression of |u| in Eq. (|120|), it can be seen that the spherical 



region of convergence defined in Eqs. (114), ( |115| ) is included in the cylindrical region of 
convergence defined in Eqs. Q125] ) and ( |126| ). 



2.6 Analytic functions of circular fourcomplex variables 

The derivative of a function f(u) of the fourcomplex variables u is defined as a function f'(u) 
having the property that 



\f(u) - /(n ) - f'(u )(u - u )\ -> as \u - u \ 0. 



(127) 



If the difference u — uq is not parallel to one of the nodal hypersurfaces, the definition in Eq. 



(128) 



(127) can also be written as 

-// \ r f( u ) ~ f( u o) 
f {u ) = hm . 

u^u u — Uq 

The derivative of the function f(u) = u m , with m an integer, is f'(u) = mu m ~ l , as can be 
seen by developing u m = [uq + (u — uq)]™ 1 as 



u ■ 



E 

p=0 



ml 



p\(m — p)\ 



u '(u-u ) p , 



(129) 



and using the definition (127). 

If the function f'(u) defined in Eq. ( |127[ ) is independent of the direction in space along 
which u is approaching uq, the function f(u) is said to be analytic, analogously to the case 
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of functions of regular complex variables. |7| The function u m , with m an integer, of the 
fourcomplex variable u is analytic, because the difference u m — u™ is always proportional to 
u — uq, as can be seen from Eq. ( |129| ). Then series of integer powers of u will also be analytic 
functions of the fourcomplex variable u, and this result holds in fact for any commutative 
algebra. 

If an analytic function is defined by a series around a certain point, for example u = 0, 

as 

oo 

/(«) = ( 13 °) 

k=0 

an expansion of f(u) around a different point uo, 

oo 

f(u) = Y,Ck(u-u ) k , (131) 

k=0 

can be obtained by substituting in Eq. (|130| ) the expression of u k according to Eq. (|129| ). 
Assuming that the series are absolutely convergent so that the order of the terms can be 
modified and ordering the terms in the resulting expression according to the increasing powers 
of u — uo yields 

00 (h 4- IV 

f( u ) = E -TUr 1 ^^ - ^o) fc - (132) 
k,i=o KA - 



Since the derivative of order k at u = uq of the function f(u) , Eq. ( 130 ), is 



/ w (uo) = g(* + 01 OJW1i i j (133) 

1=0 L 

the expansion of f(u) around u = uq, Eq. ( |132D , becomes 



00 

/(^E/'Mt"-^. ( 134 ) 

fc=o K - 

which has the same form as the series expansion of 2-dimensional complex functions. The 



relation (|134| ) shows that the coefficients in the series expansion, Eq. ( 131 ), are 



c k = ^/ (fe) K). (135) 

The rules for obtaining the derivatives and the integrals of the basic functions can be ob- 
tained from the series of definitions and, as long as these series expansions have the same form 
as the corresponding series for the 2-dimensional complex functions, the rules of derivation 
and integration remain unchanged. 
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If the fourcomplex function f(u) of the fourcomplex variable u can be expressed in terms 
of the real functions P(x, y, z, t), Q(x, y, z, t),R(x, y, z, t), S(x, y, z, t) of real variables x, y, z, t 
as 



f(u) = P(x, y, z, t) + aQ(x, y, z, t) + f3R(x, y, z, t) + -yS(x, y, z, i), (136) 

then relations of equality exist between partial derivatives of the functions P, Q, R, S. These 
relations can be obtained by writing the derivative of the function / as 



lim 



u-mo Ax + aAy + /3Az + jAt 



OP dP A dP . dP A 
—Ax + —Ay + — Az + —At 

ox oy oz at 



SdQ . dQ dQ dQ \ f dR . OR . OR . OR . N 
+a ^Ax + -^Ay + -^Az + -^At + (3 ( ^Ax + —Ay + —Az + —At 
ox oy oz at J \ ox oy oz ot , 



,'dS . d S . dS . OS . 
+7 —Ax + —Ay + — Az + —At 
\ox oy oz ot 



where the difference appearing in Eq. ( |128| ) is u — uq = Ax + aAy + (3Az + ^At. These 
relations have the same form for all systems of hypercomplex numbers discussed in this work. 

For the present system of hypercomplex numbers, the relations between the partial deriva- 
tives of the functions P, Q, R, S are obtained by setting succesively in Eq. fll37| ) Ax — ► 
0, Ay = Az = At = 0; then Ay -> 0, Ax = Az = At = 0; then Az -> 0, Ax = Ay = At = 0; 
and finally At — ► 0, Ax = Ay = Az = 0. The relations are 

dP = dQ = dR = di 

dx dy dz dt ' 

dQ = _dP = _8S = OR 

dx dy dz dt ' 

dR = _dS = _8P = dQ 
dx dy dz dt ' 

dS_ = dR = dQ = dP 
dx dy dz dt ' 



The relations ( 138[ )-( p!4l[ ) are analogous to the Riemann relations for the real and imagi- 



nary components of a complex function. It can be shown from Eqs. ( p.38[ )-( p^4l| ) that the 
component P is a solution of the equations 

Q2p d 2p ^ Q2p Q2p ^ Q2p Q2p ^ Q2p Q2p ^ 

dx 2 dy 2 ' dx 2 dz 2 ' dy 2 dt 2 ' dz 2 dt 2 

d 2p d 2 P _ {] d 2 P d 2 P _ {) 

dx 2 dt 2 ' dy 2 dz 2 
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(137) 



and the components Q,R, S are solutions of similar equations. 



As can be seen from Eqs. (142)-( |l43| ), the components P, Q, R, S of an analytic function of 
circular fourcomplex variable are harmonic with respect to the pairs of variables x, y; x, z; y, t 
and z, t, and are solutions of the wave equation with respect to the pairs of variables x, t and 
y, z. The components P, Q, R, S are also solutions of the mixed-derivative equations 

Q2p d 2p Q2p Q2p Q2p Q2p 



dxdy dzdt' dxdz dydt' dxdt dydz' 
and the components Q,R, S are solutions of similar equations. 



(144) 



2.7 Integrals of functions of circular fourcomplex variables 

The singularities of circular fourcomplex functions arise from terms of the form l/(u — uo) m , 
with m > 0. Functions containing such terms are singular not only at u = uq, but also at all 
points of the two-dimensional hyperplanes passing through uq and which are parallel to the 
nodal hyperplanes. 

The integral of a circular fourcomplex function between two points A, B along a path 
situated in a region free of singularities is independent of path, which means that the integral 
of an analytic function along a loop situated in a region free from singularities is zero, 

f(u)du = 0, (145) 

r 

where it is supposed that a surface £ spanning the closed loop T is not intersected by any 
of the two-dimensional hyperplanes associated with the singularities of the function f(u). 



Using the expression, Eq. ( |136|) , for f(u) and the fact that du = dx + ady + [5dz + jdt, the 
explicit form of the integral in Eq. (|145j ) is 

j> f{u)du = <j> [{Pdx - Qdy - Rdz + Sdt) + a{Qdx + Pdy + Sdz + Rdt) 

+(3{Rdx + Sdy + Pdz + Qdt) + j(Sdx - Rdy - Qdz + Pdt)]. (146) 

If the functions P,Q,R,S are regular on a surface S spanning the loop T, the integral along 
the loop r can be transformed with the aid of the theorem of Stokes in an integral over the 
surface E of terms of the form dP/dy + dQ/dx, dP/dz + dR/dx, dP/dt — dS/dx, dQ/dz — 
dR/dy, dQ/dt + dS/dy, dR/dt + dS/dz and of similar terms arising from the a, (3 and 7 
components, which are equal to zero by Eqs. ( |138|) -( |i~4"l|) , and this proves Eq. ( |145|) . 
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The integral of the function (u — uo) m on a closed loop T is equal to zero for m a positive 
or negative integer not equal to -1, 

<j> (u — uo) m du = 0, m integer, m ^ —1. (147) 

This is due to the fact that J(u — uo) m du = (u — uo) m+1 /{m + 1), and to the fact that the 
function (u — uo) m+1 is singlevalued for m an integer. 



The integral § du/(u — uq) can be calculated using the exponential form (|70|) 



u uq = p exp ( ol———— -\- [3——^—— + 'y hi tan ip ) , (148) 



so that 



du dp a + 3 , , a — 3 , , , , . , . 

r + —--#+ -——d X + 7d In tanV- (149) 



it — no p 2 2 
Since p and In tan -0 are singlevalued variables, it follows that §^dp/p = 0, fp din tan *0 = 0. 
On the other hand, (p and x are cyclic variables, so that they may give a contribution to 
the integral around the closed loop T. Thus, if C+ is a circle of radius r parallel to the £Ov 
plane, and the projection of the center of this circle on the £Ov plane coincides with the 
projection of the point uq on this plane, the points of the circle C+ are described according 
to Eqs. (H)-(||) by the equations 

£ = £o + r sin ip cos (p, v = v q + r sin ip sin 0, r = tq + r cos ^ cos x, 

C = Co + r cos ^ sin x, (150) 

for constant values of x and ip, ip 0, 7r/2, where «o = + «yo + /^o + 7*0 1 and Coi ^0) T 0i Co 
are calculated from xo,yo,^0)*o according to Eqs. (|l~8|). Then 

/ — = 7r(a + /3). (151) 

JC+ U — Uq 

If C_ is a circle of radius r parallel to the rOC plane, and the projection of the center of 
this circle on the tOC, plane coincides with the projection of the point uq on this plane, the 
points of the circle C_ are described by the same Eqs. (|15C| ) but for constant values of <p 
and ip, ip 7^ 0, tt/2 Then 

JC- U — Uq 
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The expression of § r du/(u — uq) can be written as a single equation with the aid of a 
functional int(M, C) defined for a point M and a closed curve C in a two-dimensional plane, 
such that 

, 1 if M is an interior point of C, 
mt(M,C) = { (153) 

if M is exterior to C. 



f du 

f =ir(a + P) mt(u t v ,Tt v ) + ir(a - (3) int(ii 0r o T r c), (154) 

Jr u — Un 



With this notation the result of the integration along a closed path T can be written as 

du 
/r u - u 

where uq^ v , uq t( ^ and T^, are respectively the projections of the point uq and of the loop 
r on the planes and tC,. 

If f(u) is an analytic circular fourcomplex function which can be expanded in a series as 



written in Eq. (131), and the expansion holds on the curve T and on a surface spanning T, 
then from Eqs. ( |147| ) and ( |154|) it follows that 

f / u)du _ + int(u ^,r f „) + (a - 0) int(ti rC> r r ()] /(u ), (155) 
Jr m - u 

where r^„,r T ^ are the projections of the curve T on the planes £v and respectively tQ, as 
shown in Fig. 2. Substituting in the right-hand side of Eq. ( |155 ) the expression of /(u) in 
terms of the real components P, Q, R, S, Eq. ( |136| ), yields 

f(u)du 
/r u — Uq 



r fWdu = + + + + + . n r 
Jr « - «o 

+tt[-(1 - 7 )(Q - R) + (a - f3)(P - S)] mt(u 0TO T T(: ), (156) 



where P, Q, R, S are the values of the components of / at u = uq. 

If f{u) can be expanded as written in Eq. ( |131| ) on T and on a surface spanning T, then 
from Eqs. (|147|) and (154) it also results that 



/ 1 f{U)d l + i =—M + P) mt(n ^,r^) + (a-/3) int(n 0rC , r rf )] f (m) (u ), (157) 
Jr [u — uq) m+i m! 



where it has been used the fact that the derivative f^ m \uo) of order m of f(u) at 



U = Uq IS 



related to the expansion coefficient in Eq. (131) according to Eq. ( |135| ). 



If a function f(u) is expanded in positive and negative powers of u — Uj, where Uj are 
circular fourcomplex constants, j being an index, the integral of / on a closed loop F is 
determined by the terms in the expansion of / which are of the form a,j/(u — Uj), 

/(„) = ... + 53 + (158) 
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Then the integral of / on a closed loop T is 

<j) f(u)du = ir(a + j3) ^ int(uj^ v ,T^ v )aj + ir(a — (3) ^ int^vo ^rc) a j- (159) 
r j j 

2.8 Factorization of circular fourcomplex polynomials 

A polynomial of degree m of the circular fourcomplex variable u = x + ay + (3z + "yt has the 
form 

P m (u) = u m + am 171 ' 1 + ■■■ + a m -iu + a m , (160) 

where the constants are in general circular fourcomplex numbers. 

It can be shown that any circular fourcomplex polynomial has a circular fourcomplex 
root, whence it follows that a polynomial of degree m can be written as a product of m 
linear factors of the form u — Uj, where the circular fourcomplex numbers Uj are the roots of 
the polynomials, although the factorization may not be unique, 

m 

Pm(u) = Y[(u- Uj ). (161) 
3=1 

The fact that any circular fourcomplex polynomial has a root can be shown by considering 
the transformation of a fourdimensional sphere with the center at the origin by the function 
u m . The points of the hypersphere of radius d are of the form written in Eq. (|80|), with d 
constant and < <fi < 2tt, < x < 2vr 7 < -0 < 7r/2. The point u m is 

u m = d m exp |ara^ + Pm^^j [cos(V - vr/4) + 7 sin(V - vr/4)] m . (162) 



It can be shown with the aid of Eq. (|84|) that 

= M, (163) 



, 4> + x . n4>-x 

uexp a — h p — - — 



so that 



[cos(V> — 7r/4) + 7sin(?/> — vr/4)] m exp ^am ^ ^ - + (3m^ 



= | [cos(^ - vr/4) + 7 sin(^ - vr/4)] m | . (164) 
The right-hand side of Eq. (|164|) is 



(cose + 7sine) m | 2 = £ C|* cos 2m - 2fc e sin 2fc e, (165) 

fe=0 
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where e = ifi — 7r/4, and since Cf^ > it can be concluded that 

|(cose + 7sine) m | 2 > 1. (166) 

Then 

d m < \u rn \ < 2 ( m - 1 )/ 2 ( r, (167) 

which shows that the image of a four-dimensional sphere via the transformation operated by 
the function u m is a finite hypersurface. 
If v! = u m , and 

v! = d'[cos(^' - vr/4) + 7 sin(^ - vr/4)] exp (c fi- | *' + Z ^'"*' ) , (168) 

then 

0' = ttk/), x = B*Xi tan?/> = tan" 1 -;/;. (169) 

Since for any values of the angles (j)',^,^' there is a set of solutions 0, x, ip of Eqs. (|169| ), 
and since the image of the hypersphere is a finite hypersurface, it follows that the image of 
the four-dimensional sphere via the function u m is also a closed hypersurface. A continuous 
hypersurface is called closed when any ray issued from the origin intersects that surface at 
least once in the finite part of the space. 

A transformation of the four-dimensional space by the polynomial P m (u) will be con- 
sidered further. By this transformation, a hypersphere of radius d having the center at the 
origin is changed into a certain finite closed surface, as discussed previously. The transfor- 
mation of the four-dimensional space by the polynomial P m {u) associates to the point u = 
the point /(0) = a m , and the image of a hypersphere of very large radius d can be repre- 
sented with good approximation by the image of that hypersphere by the function u m . The 
origin of the axes is an inner point of the latter image. If the radius of the hypersphere is 
now reduced continuously from the initial very large values to zero, the image hypersphere 
encloses initially the origin, but the image shrinks to a m when the radius approaches the 
value zero. Thus, the origin is initially inside the image hypersurface, and it lies outside the 
image hypersurface when the radius of the hypersphere tends to zero. Then since the image 
hypersurface is closed, the image surface must intersect at some stage the origin of the axes, 
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which means that there is a point u\ such that f(ui) = 0. The factorization in Eq. ( |161| ) 
can then be obtained by iterations. 

The roots of the polynomial P m can be obtained by the following method. If the constants 
in Eq. ( |160D are a; = a/o + olclh + 0ai 2 + 70/3, and with the notations of Eq. ( |123j ), the 
polynomial P m (u) can be written as 



P m = Y. 2 {m -° /2 (eiAi + eiin)(ei£ + e x i 
1=0 

m 

+ E 2 {m - /)/2 (e 2 ^ 2 + e 2 A l2 )(e 2 r + e 2 () m ~ l , (170) 
1=0 

where the constants Aik,Aik,k = 1,2, are real numbers. Each of the polynomials of degree 
m in ei£ + e\v, e 2 T + e 2 Q in Eq. ( |170| ) can always be written as a product of linear factors of 
the form ei(£ — £ p ) +e\{v — v p ) and respectively e 2 (r — r p ) + e 2 (£ — Cp)> where the constants 

£p,v p ,T p ,(p are real, 

m m 

2 (— 0/2 (eiAn + gl I a )( ei ^ + e lV ) m - 1 = J] 2 m / 2 { ei (£ - Q +et(v- v p )} , (171) 
1=0 p=i 



J] 2( m -^ 2 (e 2 ^ 2 + e 2 i Z2 )(e 2 r + e 2 C) m "' = II 2 ro/2 {e 2 (r - r p ) + e 2 (C - C P )} • (172) 
z=o p=i 



Due to the relations ( |11S| ), the polynomial P m (u) can be written as a product of factors 
of the form 

m 

P m (u) = II 2 m ' 2 { ei (£ - £ p ) + ~ei{v - v p ) + e 2 (r - r p ) + e 2 (( - ( p )} . (173) 
p=i 



This relation can be written with the aid of Eq. ( |116[ ) in the form (161), where 



Ur. 



2(ei£ p + ei^p + e 2 r p + e 2 C P ). (174) 



The roots ei£ p + eit> p and e 2 T p + e 2 Cp defined in Eqs. (|17lj ) and respectively (|172| ) may 
be ordered arbitrarily. This means that Eq. ( [L74D gives sets of m roots u\,...,u m of the 
polynomial P m (u), corresponding to the various ways in which the roots ei£ p + e\v p and 
e 2 T p + e 2 £ p are ordered according to p for each polynomial. Thus, while the hypercomplex 
components in Eqs. ( |171|) , Eqs. ( |172| ) taken separately have unique factorizations, the 
polynomial P m (u) can be written in many different ways as a product of linear factors. The 
result of the circular fourcomplex integration, Eq. ( |159|) , is however unique. 
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If, for example, P(u) = u 2 + l, the possible factorizations are P = {u — ei — e2)(u + ei + e.2) 
and P = (u — &\ + e2){u + e\ — £2), which can also be written as u 2 + 1 = (u — a)(u + a) or 
as u 2 + 1 = {u — (3)(u + (3). The result of the circular fourcomplex integration, Eq. ( |159| ), is 
however unique. It can be checked that (±ei ± e^) 2 = —e\ — &2 = — 1- 



2.9 Representation of circular fourcomplex numbers by irre- 
ducible matrices 

If T is the unitary matrix, 
/ 



V 



1 

V2 








1 





1 


1 





1 








1 





1 

V2 


1 

V2 






(175) 



(176) 



it can be shown that the matrix TUT 1 has the form 

1 f Vi 
\ F 2 

where C/ is the matrix in Eq. ( ]4l| ) used to represent the circular fourcomplex number u. In 
Eq. (176), V\, V2 are the matrices 



Vi 





(177) 



In Eq. (176), the symbols denote the matrix 




(178) 



The relations between the variables x + t, y + z,x — t,y — z for the multiplication of circular 
fourcomplex numbers have been written in Eqs. (|34|)-(|37|). The matrix TUT~ 1 provides an 
irreducible representation [|| of the circular fourcomplex number u in terms of matrices with 
real coefficients. 
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3 Hyperbolic Complex Numbers in Four Dimen- 
sions 

3.1 Operations with hyperbolic fourcomplex numbers 

A hyperbolic fourcomplex number is determined by its four components (x,y,z,t). The 
sum of the hyperbolic fourcomplex numbers (x,y,z,t) and (x', y', z' , t') is the hyperbolic 
fourcomplex number (x + x' , y + y' , z + z',t + t'). The product of the hyperbolic fourcomplex 
numbers (x, y, z, t) and (x 1 , y', z', t') is defined in this work to be the hyperbolic fourcomplex 
number (xx' + yy' + zz' + tt' , xy' + yx' + zt' + tz', xz' + zx' + yt' + ty', xt' + tx' + yz' + zy'). 

Hyperbolic fourcomplex numbers and their operations can be represented by writing the 
hyperbolic fourcomplex number (x,y, z,t) as u = x + ay + j3z + 'yt, where a, [3 and 7 are 
bases for which the multiplication rules are 

a 2 = 1, f3 2 = 1, 7 2 = 1, a/3 = f3a = 7, a 7 = 7a = 0, /3 7 = 7 /3 = a. (179) 

Two hyperbolic fourcomplex numbers u = x + oty + @z + 'yt,u' = x' + ay' + + 7^ are equal, 
u = u' , if and only if x = x',y = y',z = z',t = t'. If it = x+ay+(5z+^t, u' = x 1 '+ay' '+j3z' '+jt' 
are hyperbolic fourcomplex numbers, the sum u + u' and the product uu' defined above can 
be obtained by applying the usual algebraic rules to the sum (x + ay + (3z + jt) + (x' + ay' + 
f3z' + 'yt') and to the product (x + ay + (3z + yt)(x' + ay' + (3z' + jt'), and grouping of the 
resulting terms, 



uu = xx +yy +zz +tt' +a(xy +yx +zt' +tz')+f3{xz +zx +yt' +ty')+^(xt' +tx +yz +zy')(l%l) 
If u,u',u" are hyperbolic fourcomplex numbers, the multiplication is associative 



u + u' = x + x' + a(y + y') + [3{z + z') + j(t + t'), 



(180) 



(uu')u" = u{u'u") 



(182) 



and commutative 



uu = uu, 



(183) 
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as can be checked through direct calculation. The hyperbolic fourcomplex zero is + a ■ + 
(3 ■ + 7 • 0, denoted simply 0, and the hyperbolic fourcomplex unity is l + a- + /3-0 + 7-0, 
denoted simply 1. 

The inverse of the hyperbolic fourcomplex number u = x + ay + (3z + jt is a hyperbolic 
fourcomplex number v! = x' + ay' + f3z' + jt' having the property that 

uu' = 1. (184) 

Written on components, the condition, Eq. ( |184| ), is 

xx' + yy' + zz' + tt' = 1, 



185) 



yx' + xy' + tz' + zt' = 0, 

zx' + ty' + xz' + yt' = 0, 

tx' + zy' + yz' + xt' = 0. 

The system ( |185| ) has the solution 

, x(x 2 — y 2 — z 2 — t 2 ) + 
x = i — , (186) 

(187) 



x{x 2 - 


-y 2 - 


- z 2 - t 2 ) + 2yzt 






j 


y{-x 2 


+ y 2 


- z 2 - t 2 ) + 2xzt 






V 


z{—x 2 


-v 2 


+ z 2 -t 2 ) + 2xyt 






V 


t(-x 2 ■ 


-y 2 - 


- z 2 + t 2 ) + 2xyz 



(188) 

t'= ^ ~ y — v ZJL1^M1 } (189) 

provided that v ^ 0, where 

v = X A + y A + z 4 + t 4 - 2(x 2 y 2 + x 2 z 2 + x 2 t 2 + y 2 z 2 + y 2 t 2 + z 2 t 2 ) + 8xyzt. (190) 
The quantity v can be written as 

v = ss s s , (191) 

where 

s = x + y + z + t, s = x — y + z — t, s" = x + y — z — t, s'" = x — y — z + t. (192) 

The variables s,s',s",s"' will be called canonical hyperbolic fourcomplex variables. Then a 
hyperbolic fourcomplex number u = x + ay + (3z + yt has an inverse, unless 

8 = 0, or s' = 0, or s" = 0, or s'" = 0. (193) 
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For arbitrary values of the variables x, y, z, t, the quantity v can be positive or negative. 
If v > 0, the quantity fi = u 1 ^ will be called amplitude of the hyperbolic fourcomplex 
number x + ay + (3z + "ft. The normals of the hyperplanes in Eq. ( |193| ) are orthogonal 
to each other. Because of conditions (193) these hyperplanes will be also called the nodal 



hyperplanes. It can be shown that if uu' = then either u = 0, or u' = 0, or q, q' belong to 
pairs of orthogonal hypersurfaces as described further. Thus, divisors of zero exist if one of 
the hyperbolic fourcomplex numbers u, u' belongs to one of the nodal hyperplanes and the 
other hyperbolic fourcomplex number belongs to the straight line through the origin which 
is normal to that hyperplane, 

x + y + z + t = 0, and x' = y' = z = t' , (194) 

or 

x — y + z — t = 0, and x = —y = z = —t, (195) 

or 

x + y — z — t = 0, and x = y = —z = —t', (196) 

or 

x — y — z + t = 0, and x = —y = — z = t'. (197) 

Divisors of zero also exist if the hyperbolic fourcomplex numbers u, u' belong to different 
members of the pairs of two-dimensional hypersurfaces listed further, 

x + y = 0, z + t = and x' - y' = 0, z' - t' = 0, (198) 

or 

x + z = 0, y + 1 = and x' - z' = 0, y' - t' = 0, (199) 

or 

y + z = 0, x + t = and y' - z' = 0, x' - t' = 0. (200) 
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3.2 Geometric representation of hyperbolic fourcomplex num- 
bers 

The hyperbolic fourcomplex number x + ay + (3z + jt can be represented by the point A of 
coordinates (x, y, z, t). If O is the origin of the four-dimensional space x, y, z, t, the distance 
from A to the origin O can be taken as 

d 2 = x 2 + y 2 + z 2 + t 2 . (201) 

The distance d will be called modulus of the hyperbolic fourcomplex number x + ayA- 0z+jt, 
d = \u\. 

liu = x + ay + Pz + ^t,ui = x\ + ay\ +(3z\ +jh, u 2 = x 2 + ay 2 + /3z 2 + jt 2 , and u = u\u 2 , 
and if 

Sj — x j~\~y j~\~Zj-\~t j i Sj — Xj yj~\~Zj tj^ Sj — xj~\~yj Zj tj, Sj — xj yj zj~\~tj^j — 1, 2, (202) 
it can be shown that 

s = sis 2 , s = SiS 2 , s = SiS 2 , s = Si s 2 . (203) 

The relations (|203| ) are a consequence of the identities 

(xix 2 + yiy 2 + z 1 z 2 + tit 2 ) + (xiy 2 + y\x 2 + Z\t 2 + t x z%) 

+ (xxz 2 + z x x 2 + yit 2 + tiy 2 ) + {x 1 t 2 + t\x 2 + y\z 2 + z x y 2 ) 

= (xi + yi + zx + h)(x 2 + y 2 + z 2 + 1 2 ) (204) 

{xxx 2 + yiy 2 + z x z 2 + t\t 2 ) - (x\y 2 + y\x 2 + z\t 2 + t\z 2 ) 

+ (x 1 z 2 + z x x 2 + yit 2 + tiy 2 ) - (xih + hx 2 + y\z 2 + ziy 2 ) 

= - yi + zi- h)(x 2 -y 2 + z 2 - t 2 ) (205) 

(xtx 2 + yiy 2 + zi£2 + ht 2 ) + (xiy 2 + Vix 2 + £ii 2 + M 2 ) 

-(xiz 2 + zix 2 + yit 2 + tiy 2 ) - (xit 2 + t\x 2 + y x z 2 + ziy 2 ) 

= {x\ + y\-zi- h)(x 2 + y 2 - z 2 - t 2 ) (206) 
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(X1X2 + yiV2 + ziz 2 + ht 2 ) - (x\y 2 + y\x 2 + z x t 2 + hz 2 ) 

-(x 1 z 2 + z x x 2 + yit 2 + f 1,2/2) + (^1*2 + M2 + yi^2 + 21^/2) 

= (si - 2/1 - z\ + h)(x 2 -y 2 - z 2 + t 2 ) (207) 



A consequence of the relations (|203| ) is that if u = u\u 2 , then 

v = u x u 2 , (208) 

where 

^• = ^•444', j = l,2. (209) 

The hyperbolic fourcomplex numbers 

l + a + /3 + 7 , l-a + /3-7 „ 1 + a - (3 - 7 ,„ 1 -a-/3 + 7 
e= 5 ,e= i ,e = ,e = (210) 

are orthogonal, 

ee ' = 0, ee" = 0, ee'" = 0, e'e" = 0, e'e'" = 0, e"e"' = 0, (211) 
and have also the property that 

e 2 = e> e '2 = e / > e "2 = e « e ///2 = e /// (212) 

The hyperbolic fourcomplex number u = x + ay + (3z + 7i can be written as 

x+ay+/3z+7t = (x+y+z+t)e+(j;-y+2:-t)e / + (x+y-z-t)e" + (x-y-z+t)e / ", (213) 



or, by using Eq. ( 192 ), 

lli II II 1 W /// /01 /A 

■u = se + se+se + s e . (^14) 
The ensemble e,e',e",e"' will be called the canonical hyperbolic fourcomplex base, and Eq. 



(214) gives the canonical form of the hyperbolic fourcomplex number. Thus, if Uj = Sje + 
s'je' + s"e" + s'j'e'", j = 1,2, and u = u x u 2 , then the multiplication of the hyperbolic 
fourcomplex numbers is expressed by the relations ( j203|) . The moduli of the bases e,e',e", e'" 



arc 



M = \, W\ = \, le'1 = \, k"| = \. (215) 
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The distance d, Eq. fl20ip , is given by 

j2 1 ( 2 , /2 i 112 , ///2\ 

a=-(^s+s+s +s J. 
The relation (|216j) shows that the variables s,s',s", s'" can be written as 



s = 2d cos -0 cos cp, s' = 2d cos ip sin 0, s = 2<isin0cosx> s = 2<isin0sinx, 



(216) 



(217) 



where is the azimuthal angle in the s, s' plane, < cp < 2tt, x is the azimuthal angle in 
the s",s"' plane, < x < 27r, and ^ is the angle between the line OA and the plane ss', 
< ip < 7r/2. The variables x,y,z,t can be expressed in terms of the distance d and the 
angles 0, X> "0 as 



x = ^ (cos cos + cos sin cp + sin ^ cos x + sin sin x) > 
y = f (cos ^ cos cp — cos sin cp + sin -0 cos x — sin sin x) , 
z = | (cos -0 cos + cos ip sin — sin ip cos x — sin ip sin x) , 



(218) 



t = i (cos ip cos — cos ip sin — sin ip cos x + sin ip sin x) • 



If u = U1U2, and the hypercomplex numbers ui, U2 are described by the variables d±,cpi,xi, 0i 



and respectively c^, 02, X2> V^) then from Eq. (|217| ) it results that 



tan cp = tan cp\ tan 02 , tan x = tan xi tan X2 ; 

tan 2 ip sin 2x tan 2 ip\ sin 2xi tan 2 ip2 sin 2x2 



(219) 



sin 2cp sin 20i sin 202 

The relation ( |208] ) for the product of hyperbolic fourcomplex numbers can be demon- 
strated also by using a representation of the multiplication of the hyperbolic fourcomplex 
numbers by matrices, in which the hyperbolic fourcomplex number u = x + ay + (3z + jt is 
represented by the matrix 



(220) 



^ x y z t ^ 
y x t z 
z t x y 
\t z y x ) 

The product u = x + ay + (3z + jt of the hyperbolic fourcomplex numbers u\ = x\ + ayi + 
(3z\ + 7^1, U2 = X2 + a?/2 + (3z2 + 7*2, can be represented by the matrix multiplication 



At A 



1^2- 



(221) 
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It can be checked that the determinant det(-A) of the matrix A is 

det^ = v. (222) 

The identity ( |208| ) is then a consequence of the fact the determinant of the product of matrices 
is equal to the product of the determinants of the factor matrices. 

3.3 Exponential form of a hyperbolic fourcomplex number 

The exponential function of a fourcomplex variable u and the addition theorem for the 
exponential function have been written in Eqs. (Hj) and ([II]). If u = x + ay + 0z + jt, 
then exp u can be calculated as exp u = exp x ■ exp(ay) • exp(/3z) • exp(7i) . According to Eqs. 

(ED, 



a 2m = 1, a 2m+1 = a, (3 2m = 1, f3 2m+1 = (3, 7 2m = 1, 7 2m+1 = 7, (223) 
where m is a natural number, so that exp(ay), exp(/3z) and exp(jt) can be written as 

exp(ay) = coshy+asinhy, exp(/?z) = cosh y+/3sinhz, exp ( 7 t) = cosh t+jsinht. (224) 



From Eqs. (224) it can be inferred that 

(cosh t + a sinh t) m = cosh mt + a sinh mt , (cosh t + (3 sinh t) m = cosh mi + /3 sinh mt, 
(cosh t + 7 sinh t) m = cosh mt + 7 sinh mt. (225) 

The hyperbolic fourcomplex numbers u = x + ay + /3z + jt for which s = x + y + z + t> 
0, s' = x — y + z — t > 0, s" = x + y — z — t > 0, s" = x — y — z + t > can be written in 
the form 

x + ay + (3z + jt = e xl+ayx+ P zl+n < tl . (226) 

The conditions s = x + y + z + t>0, s' = x — y + z — t > 0, s" = x + y — z — i > 0, s'" = 
x — y — z + t > correspond in Eq. (|2 1 7| ) to a range of angles < < 7r/2,0 < x < 
7r/2,0 < if) < ir/2. The expressions of X\, yi, z%, t\ as functions of x,y,z,t can be obtained 
by developing e ayi , e@ Zl and e 7 * 1 with the aid of Eqs. ( |224f ) , by multiplying these expressions 
and separating the hypercomplex components, 

x = e xi (cosh yi cosh z\ cosh t\ + sinh yi sinh 21 sinh ii ) , (227) 
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y = e ai (sinhyi coshzi coshii + coshyi sinh^i sinhti), 



(228) 



z = e Xl (cosh yi sinh z\ cosh t\ + sinh y\ cosh zi sinh t\) , (229) 
i = e Xl (sinh yi sinh zi cosh ii + cosh yi cosh #i sinh t\) , (230) 

It can be shown from Eqs. (|227|)-(|230|) that 

1 1 ss" 1 ss' 1 ss'" 

xi = — hi(ss s's"), y± = — In , „, , zi = — In , t\ = — In -r-^- (231) 

The exponential form of the hyperbolic fourcomplex number u can be written as 
1 , ss" 1 ss' 1 , ss'" 



u 



^(r'^ + Z'^ + ^ln— j, (232) 
where 

H = {ss's's'") 1 ^. (233) 

The exponential form of the hyperbolic fourcomplex number u can be written with the aid 
of the relations (|217|) as 



(\ 1 1 sm2<^> 1 tanv\ , .» 

it = a exp -a In ■ h -p In — - — h -7 In . (234) 

\4 tancptanx 4 tan z ^ sm 1\ 4 tan cp J 

The amplitude /i can be expressed in terms of the distance d with the aid of Eqs. ( |217| ) as 

H = d sin 1/2 2ip sin 1/4 2(f> sin 1/4 2 X . (235) 

The hypercomplex number can be written as 

u = dsin 1/2 2V>sin 1/4 2</)sin 1/4 2 X exp (\a\n p + i/?ln i^-^ + ^7 In ,(236) 

V 4 tan cp tan x 4 tan z ip sm 2% 4 tan q> ) 

which is the trigonometric form of the hypercomplex number u. 



3.4 Elementary functions of a hyperbolic fourcomplex vari- 
able 

The logarithm u\ of the hyperbolic fourcomplex number u, u\ = In it, can be defined for 
s > 0, s' > 0, s" > 0, s'" > as the solution of the equation 

u = e Ul , (237) 
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for u\ as a function of u. From Eq. ( |232| ) it results that 
Inu 



1. 1 . ss" 1 , ss' 1 , ss w 

In At + +-aln —— + -/31n — — + -7 In — . 



4 r 4 sV" 4' s'V" 4 s's" 
Using the expression in Eq. (|234| ) , the logarithm can be written as 

1 



1, 1 . 1 1 „ . sm 2</> 1 tan x 
- In /i + -ram ■ h -p in s— - — : — h -7m -. 

4 4 tan tan % 4 tan z sm 2% 4 tan <p 



In it 



It can be inferred from Eqs. (f238|) and ( |203|) that 



(238) 



(239) 



ln(iti«2) = lnui +lnit2. 



(240) 



The explicit form of Eq. ( 238 ) is 



ln(x + ay + (3z + jt) = -(1 + a + + 7) ln(x + y + 2 + t) 

+-(1 - a + /3 - 7) ln(x - y + 2 - i) + -(1 + a - (3 - 7) ln(x + y - z - t) 
+-(1 - a - /3 + 7) ln(x - y - z + i). 



(241) 



The relation ( 241 ) can be written with the aid of Eq. fl210|) as 



In u = e In s + e' In s' + e" In s" + e'" In s'" 



(242) 



The power function u n can be defined for s > 0, s' > 0, s" > 0, s'" > and real values of 



n as 



(243) 



It can be inferred from Eqs. (243) and Q240| ) that 



(u lU2 ) n = u\ u n 2 . 



(244) 



Using the expression (241) for Inn and the relations fl211| ) and (212) it can be shown that 



(x + ay + 0z + 7t) n = i(l + a + + j)(x + y + z + t) n + i(l - a + /? - 7)(x - y + z - t)' 
+~(1 + a- /3-7)(x + y- z-t) n + i(l-a-/3 + 7 )(x-y-z + t) n . (245) 



For integer n, the relation fl245|) is valid for any x,y,z,t. The relation (245) for n = —1 is 

1 l/l + a + /3 + 7 l-a + /3-7 1 + a - (3 - 7 l-a-/? + 7 N 



+ 



+ 



+ 



x + ay + 0z + jt A\x + y + z + t x — y + z — t x + y — z — t x — y — z + t 



.(246) 
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The trigonometric functions of the fourcomplex variable u and the addition theorems for 
these functions have been written in Eqs. (|90"|)-(|9"5|). The cosine and sine functions of the 
hypercomplex variables ay, (3z and jt can be expressed as 

cos ay = cosy, sin ay = asiny, (247) 

cos f3y = cosy, sin/?y = fls'my, (248) 

cos7y = cosy, shi7y = 7siny. (249) 

The cosine and sine functions of a hyperbolic fourcomplex number x + ay + f3z + jt can then 
be expressed in terms of elementary functions with the aid of the addition theorems Eqs. 



H), (H) and of the expressions in Eqs. (P47|)-(P4SD 



The hyperbolic functions of the fourcomplex variable u and the addition theorems for 
these functions have been written in Eqs. (|97|)-( |l00|) . The hyperbolic cosine and sine func- 
tions of the hypercomplex variables ay, (5z and "ft can be expressed as 

cosh ay = coshy, sinh ay = asinhy, (250) 

cosh/3y = coshy, sinh/?y = /3sinhy, (251) 

cosily = coshy, sinh7y = 7 sinh y. (252) 

The hyperbolic cosine and sine functions of a hyperbolic fourcomplex number x+ay+f3z+jt 
can then be expressed in terms of elementary functions with the aid of the addition theorems 
Eqs. (PG| ), ( |100D and of the expressions in Eqs. ( j250|) -( p52|) , 

3.5 Power series of hyperbolic fourcomplex variables 

A hyperbolic fourcomplex series is an infinite sum of the form 

a + ai+a 2 -\ h «H , (253) 

where the coefficients ai are hyperbolic fourcomplex numbers. The convergence of the series 



( 253 ) can be defined in terms of the convergence of its 4 real components. The convergence 



of a hyperbolic fourcomplex series can however be studied using hyperbolic fourcomplex 
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variables. The main criterion for absolute convergence remains the comparison theorem, but 
this requires a number of inequalities which will be discussed further. 

The modulus of a hyperbolic fourcomplex number u = x + ay + (3z + ^ft can be defined as 



(* 2 + y 2 + z 2 +i 2 ) 1/2 , 



(254) 



so that according to Eq. (|201[) d = \u\. Since |x| < |u|, \y\ < \u\, \z\ < \u\, \ t\ < \u\, a property 
of absolute convergence established via a comparison theorem based on the modulus of the 



series ( 253 ) will ensure the absolute convergence of each real component of that series. 

The modulus of the sum u\ + ui of the hyperbolic fourcomplex numbers u±,U2 fulfils the 
inequality 



\ui\ — \U2\\ < \ui + U2I < Wl\ + \U2\ 



For the product the relation is 



IU1W2I < 2|ui||u2| 



(255) 



(256) 



which replaces the relation of equality extant for regular complex numbers. The equality 
in Eq. ( |256| ) takes place for x\ = y\ = z\ = t\ and x%jx\ = 2/2/2/1 = zi/z\ = ^Ai- In 
particular 



\u 2 \ < 2(x 2 + y 2 + z 2 + t 2 ). 



The inequality in Eq. (|256j) implies that 



\u l \<2 l - l \u\ l . 



From Eqs. ( 256 ) and fl25q ) it results that 



\au l \ < 2 l \a\\u\ l . 



A power series of the hyperbolic fourcomplex variable u is a series of the form 



ao + a\u + a^u + • • • + a[U + 



Since 



<E 2 'i a «iki', 
1=0 



(257) 



(258) 



(259) 



(260) 



(261) 
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a sufficient condition for the absolute convergence of this series is that 

lim ^±iM < i. (262) 

Thus the series is absolutely convergent for 

M < c , (263) 

where 

c = lim (264) 



The convergence of the series ( |260| ) can be also studied with the aid of the formula (245) 
which, for integer values of I, is valid for any x, y, z, t. If ai = ai x + aai y + f3a\ z + 70/4, and 

Ai = ai x + ai y + ai z + au, (265) 

A'i = a\ x - ai y + a\ z - au, (266) 

A" = a\ x + ai y - au - au, (267) 

A'" = a\ x - a\y - a\ z + au, (268) 



it can be shown with the aid of relations ( |211| ) and (212) that 

A I A> I II A'I II HI AlH HI /lcn\ 

a/e = Aie, aie = A t e , aie = A l e , a/e = A x e , (269) 
so that the expression of the series (|260|) becomes 



J2 (A lS l e + A' lS ' l e' + A'{ s" l e" + AfV'V") , (270) 
1=0 

where the quantities s,s' ,s" , s'" have been defined in Eq. ( |192j ). The sufficient conditions for 
the absolute convergence of the series in Eq. (|270|) are that 



\Ai +1 \\s\ \A' J+1 \\s'\ \A','A\s"\ \A'/'\\s"'\ , 

lim 1 +1 1 1 < 1, lim 1 + 1 ' 1 < 1, lim 1 1 < 1, lim 1 1 < 1, (271) 

l— >oo \Ai\ l— >oo \Ai\ l^oo \A^ I l^oo \A^ I 

Thus the series in Eq. (|270|) is absolutely convergent for 



\x + y + z + t\ < c, \x — y + z — t\ < c, \x + y — z — t\ < c , \x — y — z + t\ < c" , (272) 

where 

14,1 \A'\ \A"\ \A'"\ 

c = lim d = lim -^L, c" = lim J^L, d" = lim ( 27 3) 
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The relations ( p72| ) show that the region of convergence of the series ( |270 ) is a four-dimensional 



parallelepiped. It can be shown that cq = (l/2)min(c, c', c", c'"), where min designates the 
smallest of the numbers c,c',c",d". Using Eq. fl216|), it can be seen that the circular region 



of convergence defined in Eqs. (263), (264) is included in the parallelogram defined in Eqs. 
p72D and (g7g). 



3.6 Analytic functions of hyperbolic fourcomplex variables 



The fourcomplex function f(u) of the fourcomplex variable u has been expressed in Eq. ( |136| ) 
in terms of the real functions P(x, y, z, t), Q(x, y, z, t),R(x, y, z, t), S(x, y, z, t) of real variables 
x, y, z, t. The relations between the partial derivatives of the functions P, Q, R, S are obtained 
by setting succesively in Eq. ( |137| ) Ax —* 0, Ay = Az = At = 0; then Ay — > 0, Ax = Az = 
At = 0; then Az -> 0, Ax = Ay = At = 0; and finally At -> 0, Ax = Ay = Az = 0. The 
relations are 

PlD FlU P.Q 

(274) 
(275) 
(276) 
(277) 

The relations ( p74|) -( ^77D are analogous to the Riemann relations for the real and imag- 
inary components of a complex function. It can be shown from Eqs. fl274D -( f277| ) that the 
component P is a solution of the equations 

d 2 P d 2 P _ cPP_ &P_ _ (PP_ cPP_ _ (PP_ _ 

dx 2 dy 2 ' dx 2 dz 2 ' dy 2 dt 2 ' dz 2 dt 2 

dx 2 dt 2 ' dy 2 dz 2 

and the components Q, R, S are solutions of similar equations. As can be seen from Eqs. 
( 278[ )-( p79|) , the components P,Q,R,S of an analytic function of hyperbolic fourcomplex 



dP 


_ dQ _ 


dR 


dS 


dx 


dy 


dz 


~ ~dt 


dQ 


dP 
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dR 
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dz 


~ ~dt 


OR 
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dQ 
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variable are solutions of the wave equation with respect to pairs of the variables x, y, z, t. 
The component P is also a solution of the mixed-derivative equations 

g2p ^ Q2p Q2p ^ Q2p Q2p ^ Q2p 

dxdy dzdV dxdz dydV dxdt dydz' 
and the components Q, R, S are solutions of similar equations. 



3.7 Integrals of functions of hyperbolic fourcomplex variables 

The singularities of hyperbolic fourcomplex functions arise from terms of the form l/(u — 
uo) m , with m > 0. Functions containing such terms are singular not only at u = uq, but also 
at all points of the two-dimensional hyperplanes passing through uq and which are parallel 
to the nodal hyperplanes. 

The integral of a hyperbolic fourcomplex function between two points A, B along a path 
situated in a region free of singularities is independent of path, which means that the integral 
of an analytic function along a loop situated in a region free from singularities is zero, 

f(u)du = 0, (281) 

r 

where it is supposed that a surface £ spanning the closed loop T is not intersected by any 
of the two-dimensional hyperplanes associated with the singularities of the function f(u). 



Using the expression, Eq. ( |136|) , for f(u) and the fact that du = dx + ady + fldz + 'jdt, the 
explicit form of the integral in Eq. (|281] ) is 

j> f(u)du = j> [{Pdx + Qdy + Rdz + Sdt) + a{Qdx + Pdy + Sdz + Rdt) 

+(3{Rdx + Sdy + Pdz + Qdt) + j(Sdx + Rdy + Qdz + Pdt)] . (282) 

If the functions P,Q,R,S are regular on a surface X spanning the loop T, the integral along 
the loop r can be transformed with the aid of the theorem of Stokes in an integral over the 
surface S of terms of the form dP/dy — dQ/dx, dP/dz — dR/dx, dP/dt — dS/dx, dQ/dz — 
dR/dy, dQ/dt — dS/dy, dR/dt — dS/dz and of similar terms arising from the a, (5 and 7 
components, which are equal to zero by Eqs. ( |274 )-( p77|) , and this proves Eq. fl281|) . 



The exponential form of the hyperbolic fourcomplex numbers, Eq. fl234|) , contains no 
cyclic variable, and therefore the concept of residue is not applicable to the hyperbolic four- 



complex numbers defined in Eqs. ( 179 ). 
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3.8 Factorization of hyperbolic fourcomplex polynomials 



A polynomial of degree m of the hyperbolic fourcomplex variable u = x + ay + (3z + <yt has 
the form 



P m (u) = u m + a x u m H V cim-xu + a r , 



(283) 



where the constants are in general hyperbolic fourcomplex numbers. If a m = a mx + aa my + 
Pamz + ia-mt, and with the notations of Eqs. (|192| ) and (|265[) -( p68|) applied for I = 0, 1, • • • , m 
, the polynomial P m (u) can be written as 

s m + ^is m_1 H + A m _is + A m l e + 



Prr 



+ 



, All iim-1 j , 4// // , 4// 



lllm + A 'n s lllm-1 + . . . + J 4'^_ lS '" + ,4; 



e" + 



m 

/// „/// 



Each of the polynomials of degree m with real coefficients in Eq. fl284|) can be written as a 
product of linear or quadratic factors with real coefficients, or as a product of linear factors 
which, if imaginary, appear always in complex conjugate pairs. Using the latter form for the 
simplicity of notations, the polynomial P m can be written as 

m m m m 

P m = U( s ~ s l> + IK*' - ^ + IK S " - S "V + IK 8 '" - (285) 
1=1 1=1 1=1 1=1 

where the quantities S[ appear always in complex conjugate pairs, and the same is true for 
the quantities s 1 ,, for the quantitites sf, and for the quantities s 1 ,". Due to the properties in 



Eqs. (211) and (212), the polynomial P m {u) can be written as a product of factors of the 
form 



Pm(u) = II [(« - Sl )e + (s> - s[)e' + (s" - s'l)e" + {s'" - s'l')e"'] 
1=1 

These relations can be written with the aid of Eq. ( [21 3D as 

P m {u) 

where 



Y[(u-u p ), 

p=i 



. / / 1 //„// j_ Jtt Jtf 



(286) 



(287) 



(288) 



The roots s p , s', s" s"' of the corresponding polynomials in Eq. (285) may be ordered arbi- 



-"pi "pi "p 

trarily. This means that Eq. (|288| ) gives sets of m roots u±, ...,u m of the polynomial P m (u) 



(284) 
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corresponding to the various ways in which the roots s p , sL, s'L s 1 '' are ordered according to 



p in each group. Thus, while the hypercomplex components in Eq. (284) taken separately 
have unique factorizations, the polynomial P m (u) can be written in many different ways as 
a product of linear factors. 



If P(u) 



1, the factorization in Eq. (|287|) is u 2 — 1 = (u — u\)[u — U2), where 



ui = ±e ± e' ± e" ± e'" , u 2 = — Ui, so that there are 8 distinct factorizations of u 2 — 1, 



u 2 
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= (u- 
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3.9 Representation of hyperbolic fourcomplex numbers by ir- 
reducible matrices 



If T is the unitary matrix, 



/ 1 



1 

2 

1 

2 

1 

V 2 



it can be shown that the matrix TUT 1 has the form 



TUT' 1 



x+y+z+t 






x 





y + z 







t 

x + y - z 








x — y — z + t 



(290) 



(291) 
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where U is the matrix in Eq. ( |22C| ) used to represent the hyperbolic fourcomplex number u. 
The relations between the variables x + y + z + t,x — y + z — t,x + y — z — t, x — y — z + t for 
the multiplication of hyperbolic fourcomplex numbers have been written in Eqs. ( |204| )- ( pQ7| ) . 
The matrix TUT -1 provides an irreducible representation Q of the hyperbolic fourcomplex 
number u in terms of matrices with real coefficients. 

4 Planar Complex Numbers in Four Dimensions 
4.1 Operations with planar fourcomplex numbers 

A planar fourcomplex number is determined by its four components (x,y,z,t). The sum of 
the planar fourcomplex numbers (x, y, z, t) and (a/, y' , z' ,t') is the planar fourcomplex number 
(x + x' , y + y', z + z',t + t'). The product of the planar fourcomplex numbers (x, y, z, t) and 
(x' , y' , z' ,t') is defined in this work to be the planar fourcomplex number (xx' — yt' — zz' — 
ty' , xy' + yx' — zt' — tz' , xz' + yy' + zx' — tt' ', xt' + yz' + zy' + tx'). 

Planar fourcomplex numbers and their operations can be represented by writing the 
planar fourcomplex number (x, y,z,t) asu = x + ay + j3z + yt, where a, (3 and 7 are bases 
for which the multiplication rules are 

a 2 = /3, p 2 = -1, 7 2 = -13, a(3 = f3a = 7, 07 = 7a = -1, /3y = yf3= -a. (292) 

Two planar fourcomplex numbers u = x + ay + [3z + yt, v! = x' + ay' + (3z' + yt' are equal, 
u = u', if and only if x = x', y = y' , z = z' , t = t' . If u = x+ay+(3z+yt, v! = x' +ay' +(3z' +yt' 
are planar fourcomplex numbers, the sum u + u' and the product uu' defined above can be 
obtained by applying the usual algebraic rules to the sum (x + ay + 0z + yt) + (x' + ay' + 
(3z' + yt') and to the product (x + ay + (3z + yt){x' + ay' + f3z' + yt'), and grouping of the 
resulting terms, 

u + u ' = x + x' + a(y + y') + f3{z + z) + 7 (t + t'), (293) 

uu' = xx' — yt' — zz — ty' + a(xy + yx' — zt! — tz) + (3{xz + yy' + zx' — tt') 

+y{xt' + yz' + zy +tx). (294) 
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If u,u',u" are planar fourcomplex numbers, the multiplication is associative 

(uu)u' = u(u'u") (295) 
and commutative 

uu = uu, (296) 

as can be checked through direct calculation. The planar fourcomplex zero is + a ■ + (3 ■ 
+ 7-0, denoted simply 0, and the planar fourcomplex unity is l + a- + /3-0 + 7- 0, denoted 
simply 1. 

The inverse of the planar fourcomplex number u = x+ay-\-f3z+ r yt is a planar fourcomplex 
number u' = x' + ay' + f3z' + jt' having the property that 

uu = 1. (297) 

Written on components, the condition, Eq. ( |297j ), is 

xx' — ty' — zz' — yt' = 1, 



(298) 



yx' + xy' — tz' — zt' = 0, 

zx' + yy' + xz' - tt' = 0, 

tx' + zy' + yz' + xt' = 0. 

The system fl298|) has the solution 

x , = x(x 2 + z 2 ) - z(y 2 - t 2 ) + 2xyt 
P 4 

y(x 2 - z 2 ) + t(y 2 + t 2 ) + 2xzt 
V = —± , (300) 

z , = -z(x 2 + z 2 ) + x(y 2 - t 2 ) + 2zyt 
2 P 4 
t , = t(x 2 - z 2 ) + y(y 2 + t 2 ) - 2xyz 
P 4 

provided that p ^ 0, where 

p 4 = x 4 + z 4 + y 4 + i 4 + 2x 2 z 2 + 2y 2 t 2 + 4x 2 yt - Axy 2 z + Axzt 2 - 4yz 2 t. (303) 

The quantity p will be called amplitude of the planar fourcomplex number x + ay + (3z + jt. 
Since 

P 4 = P 2 + P-, (304) 
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where 



2 / , y-t\ 2 , ( , y + t\ 2 2 (_ y-t\ 2 , (. y + t^ 2 



a planar fourcomplex number u = x + ay + (3z + 7* has an inverse, unless 

y — t y + t , 

x+^=- = 0, z + ^=- = 0, (306) 



or 



y~ l n y + t n Con^ 

x-- 7 =- = 0,z-- 7 =- = 0. (307) 



Because of conditions ( |306| )- (|307l ) these 2-dimensional hypersurfaces will be called nodal 
hyperplanes. It can be shown that if uv! = then either u = 0, or v! = 0, or one of the 
planar fourcomplex numbers is of the form x + a{x + z)j\pl + (3z — ^{x — z)/y/2 and the 
other of the form x' - a(x' + z') /y/2 + (3z' + <y(x' - z')/V2. 

4.2 Geometric representation of planar fourcomplex numbers 

The planar fourcomplex number x + ay + f3z + yt can be represented by the point A of 
coordinates (x, y, z, t). If O is the origin of the four-dimensional space x, y, z, t, the distance 
from A to the origin O can be taken as 

d 2 = x 2 + y 2 + z 2 + t 2 . (308) 

The distance d will be called modulus of the planar fourcomplex number x + ay + j3z + jt, 
d = \u\. The orientation in the four-dimensional space of the line OA can be specified with 
the aid of three angles (f), x, ij) defined with respect to the rotated system of axes 
x y—t x y—t z y+t z y+t 

The variables £, v, r, £ will be called canonical planar fourcomplex variables. The use of the 
rotated axes £, v, r, £ for the definition of the angles eft, x, ift is convenient for the expression 
of the planar fourcomplex numbers in exponential and trigonometric forms, as it will be 
discussed further. The angle <j) is the angle between the projection of A in the plane £, v 
and the 0£ axis, < <p < 2tt, \ is the angle between the projection of A in the plane r, £ 
and the Or axis, < % < 2vr, and ij) is the angle between the line OA and the plane tO(, 
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< ip < 7r/2, as shown in Fig. 1. The definition of the variables in this section is different 
from the definition used for the circular fourcomplex numbers, because the definition of the 
rotated axes in Eq. ( |309| ) is different from the definition of the rotated circular axes, Eq. 
(JT8|) . The angles <p and x wm be called azimuthal angles, the angle ip will be called planar 
angle . The fact that < ip < ir/2 means that ip has the same sign on both faces of the 
two-dimensional hyperplane vOC,. The components of the point A in terms of the distance 
d and the angles <p, x, ip are thus 

-^= + - — - = dcos <ps\n.ip, (310) 
V 2 2 

^ - = dcosxcosip, (311) 

— p + - ~t = rfsin^sin-^, (312) 
v2 2 

~~ "71 + ^1T~ = dsmxcos ^- ( 313 ) 

It can be checked that p + = y/2dsmtp, p_ = v^dcos^. The coordinates x,y,z,t in terms of 
the variables d, 0, x, ^ are 

x = — p (cos sin + cos x cos V 7 ), (314) 
V 2 

y = —i=[sm((j) + vr/4) sin ip + sin(x — tt/4) cos (315) 
V2 

2 = -^= (sin sin ?/> — sin x cos ip) , (316) 
v2 

t = —i=[— cos((p + vr/4) sin ?/> + cos(x — vr/4) cos (317) 
V2 

The angles (f>, x, ^ can be expressed in terms of the coordinates x, y, z, t as 
. z + (y + t)/V2 x+(y-t)/y/2 

sm0= , cos<p = , (318) 

P+ P+ 

-z + (y + t)/V2 x-(y- t)/y/2 

smx = , cosx = , (319) 

P- P- 

tanip = p + /p- (320) 
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The nodal hyperplanes are £Ov , for which r = 0, £ = 0, and rOC, for which £ = 0, v = 0. 
For points in the nodal hyperplane £Ov the planar angle is ip = vr/2, for points in the nodal 
hyperplane tOQ the planar angle is tp = 0. 

It can be shown that if u\ = x\ + ayi + + 7*1,^2 = ^2 + ^2/2 + ^2 + 7*2 are planar 
fourcomplex numbers of amplitudes and angles p\ , 4>\ , xi , ipi an d respectively p 2 ,(p2,X2,'>p2, 
then the amplitude p and the angles <f>, x> ^ °f the product planar fourcomplex number ui« 2 
are 



P = P1P2, 



(321) 



h +^2, X = Xi + X2, tam/> = tam/>i tamp 2 . 



(322) 



The relations (|32l|) - (|322j) are consequences of the definitions (|303| )- (]305|) , (|31S| )- (|32C| ) and of 
the identities 



(xix 2 - z 1 z 2 - yit 2 - hy 2 ) + 



(xiy 2 + yix 2 - Zit 2 - tiz 2 ) - (x±t 2 + t\x 2 + z x y 2 + y±z 2 ) 



(x 1 z 2 + zix 2 + yiy 2 - ht 2 ) + 



V2 

(xiy 2 + yix 2 - z\t 2 - t\z 2 ) + {x\t 2 + t\x 2 + ziy 2 + y\z 2 ) 



V2 



Xl + 



yi - h 

V2 



+ Ui + 



3/1 + *i 
V2 



x 2 + 



Z/2 ~*2 

72 



+ U2 + 



2/2 + *2 

V2 



(323) 



(xix 2 - z x z 2 - yit 2 - tiy 2 ) 



(xiy 2 + yix 2 - z\t 2 - t\z 2 ) - {x\t 2 + t\x 2 + ziy 2 + y\Z 2 ) 



(x 1 z 2 + ^ix 2 + yiy 2 - ht 2 ) 



V2 

{x\y 2 + y\x 2 - z x t 2 - t x z 2 ) + (x\t 2 + hx 2 + ziy 2 + yiz 2 ) 



Vi-h\ , ( yi + h 

Xl ■=— + [Zl 



V2 



V2 



V2 

y2 - 1 2 \ 2 f y2 + t 2 
x 2 ■=- + \z 2 



V2 



V2 



(324) 



(xix 2 - z\z 2 - y x t 2 - hy 2 ) + 



(xiy 2 + 2/1X2 - zit 2 - t\z 2 ) - {x\t 2 + tix 2 + z\y 2 + y\Z 2 ) 



. yi - 1 A ( , V2-h 

Xl H — ) x 2 + 



V2 



V2 



V2 

. yi + h \ ( v2 + t 2 

zi H 7^- U2 + 



V2 



V2 



(325) 



(xiz 2 + 2:1x2 + yiy2 - ht 2 ) + 



(xiy 2 + y\x 2 - z x t 2 - t\z 2 ) + {x\t 2 + tix 2 + z x y 2 + yiz 2 ) 



V2 



. yi+h\ { .y2-t 2 \ ( ,yi-t x \( .y2 + h 

zi H x 2 H + xi H — U 2 + 



V2 



V2 



V2 



(326) 
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(xix 2 - ziz 2 - y\t 2 - hy 2 ) 



(xiy 2 + yix 2 - zit 2 - tiz 2 ) - {xit 2 + t x x 2 + ziy 2 + y x z 2 ) 



V2 



yi - h \ ( yz-h 

Xl 7=— ) \X 2 - 



V2 



-(x 1 z 2 + zix 2 + yxy 2 - tit 2 ) + 



V2 J \' Z1 1 ^ A " V2 

{x\y 2 + y\x 2 - z\t 2 - t\z 2 ) + {x\t 2 + t\x 2 + z x y 2 + y 1 z 2 ) 



yi +h 



-z 2 + 



2/2 + h 



(327) 



V2 



The identities (|23|) and d324| ) can also be written as 



-z 2 + 



2/2 + 1 2 
V2 



P+ = P1+P2+, 



P- = P1-P2-, 



where 



, Vj-tA 2 , ( , Vj + t A 2 2 

^•i + nr I + U-f + ^ J I , Pl- 



(328) 

(329) 
(330) 



V -^\ +(z- V A±\ ,(331) 



) v J 72 ; \r 3 V2 ) v 72 

for j = 1,2. 

The fact that the amplitude of the product is equal to the product of the amplitudes, as 
written in Eq. ( |321| ), can be demonstrated also by using a representation of the multiplication 
of the planar fourcomplex numbers by matrices, in which the planar fourcomplex number 
u = x + ay + f3z + "ft is represented by the matrix 

/ x y z t ^ 
-t x y z 
-z —t x y 
\-y -z -t x J 

The product u = x + ay + f3z + jt of the planar fourcomplex numbers u\ = x\ + ayi + + 
7*1, u 2 = x 2 + ay 2 + (3z 2 + jt 2 , can be represented by the matrix multiplication 



A 



(332) 



A = AjA 2 . 

It can be checked that the determinant det(^4) of the matrix A is 
detA = p 4 . 



(333) 



(334) 



The identity ( |321| ) is then a consequence of the fact the determinant of the product of matrices 
is equal to the product of the determinants of the factor matrices. 
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,8m _ 


l,a 8m+1 = 


a , a 8m+2 = ^ a 8m+3 = % 


,8m+4 


= -l,a 8m - 


"5 = _ aj a 8m+6 _ a 8m+7 _ _^ 


|4m 


l,p 4m + 1 = 


P ) pton+2 = -l,^m+3 = _^ 


,8m 




7, 7 8m + 2 = -/?, 7 8m+3 = a, 


,8m+4 


_ ^,8m4 


-5 = _ 7)7 8m+6 = p n 8m+7 = 



4.3 The planar fourdimensional cosexponential functions 

The exponential function of a fourcomplex variable u and the addition theorem for the 
exponential function have been written in Eqs. ( f44[) and (fig). If -u = x + ay + /3z + 7t, then 
exp u can be calculated as exp u = exp x ■ exp(ay) • exp(/3z) ■ ex.p(jt). According to Eq. fl292| ), 

a 8 
a 

(335) 

m _ 8m+l _ 7 ^8m+2 _ _^ ,y8m+3 _ a ^ 

where n is a natural number, so that exp(ay), exp(/3z) and exp(7i) can be written as 

exp((3z) = cos z + j3 sin z, (336) 

and 

exp(ay) = f A0 {y) + af 41 (y) + (3f A2 {y) + 7/43(2/), (337) 
exp( 7 t) = / 40 (t) + 7/41 (t) - AfeW + a/ 43 (t), (338) 

where the four-dimensional cosexponential functions /40, /41, /42, /43 are defined by the series 
/ 40 (x) = l-x 4 /4! + x 8 /8!----, (339) 

/ 41 ( x ) = x - x 5 /5! + x 9 /9\ , (340) 

U 2 {x) = x 2 /2\ - x 6 /6! + x 10 /10! , (341) 

f i3 (x) = x 3 /3\ - x 7 /7l + x n /lV. . (342) 

The functions /40, /42 are even, the functions /41, ^43 are odd, 

Uo(-u) = Uo(u), /42(-«) = /42(«), /4l(-«) = -Ul(u), f4s(-u) = -f4s(u). (343) 

Addition theorems for the four-dimensional cosexponential functions can be obtained from 
the relation exp a(x + y) = exp ax ■ exp ay, by substituting the expression of the exponentials 
as given in Eq. (|337|) , 

ho{x + y) = fao{x)fao{y) - fai(x)U 3 (y) - /42 (20/42(1/) - fa3{x)fai(y), (344) 
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Ui{x + y) = f A o(x)f A1 (y) + fa(x)f m (y) - f A2 (x)f A3 (y) - h3(x)f A2 (y), (345) 
/42O + y) = fm{x)f A2 {y) + fai(x)fai(y) + /42(«)/4o(y) - f43(x)f A3 (y), (346) 
/ 43 (x + y) = /4o(x)/ 43 (y) + fa(x)fa(v) + f42{x)My) + /43(x)/ 40 (y). (347) 



For x = y the relations ( |344| )-( |347|) take the form 

/ 40 (2x) = /| (x) - /| 2 (x) - 2/ 41 (x)/ 43 (x), (348) 

f A1 (2x) = 2f 40 (x)f 41 (x) - 2f A2 (x)f A3 (x), (349) 

/ 42 (2x) = fl(x) - fUx) + 2f 40 (x)f 42 (x), (350) 

/ 43 (2x) = 2f A0 (x)f A3 (x) + 2f A1 (x)f A2 (x). (351) 
For x = -y the relations (^4^) - (|347|) and (|343|) yield 

/lo(^) " + 2/ 4 i(x)/ 43 (x) = 1, (352) 

fli(x) - fUx) - 2f A0 (x)f A2 (x) = 0. (353) 



From Eqs. (|336[ )- ( p38|) it can be shown that, for m integer, 



(cos z + sin z) m = cos mz + (3 sin mz, (354) 

and 

[Uo(y)+af4i(y)+f3fa 2 (y)+'7fa3(y)} m = fao{rny)+af A1 (my)+(3f A2 (my)+<yf A3 (my), (355) 
[/ 40 it) + 7/41 (t ) - /?/ 42 (t) + «/ 43 (t )] m = f A0 (mt) + 7 / 41 (mt) - /3/ 42 (mt) + af A3 (mt) . (356) 
Since 

(a - 7 ) 2m = 2 m , (a - 7 ) 2m+1 = 2 rn (a - 7), (357) 
it can be shown from the definition of the exponential function, Eq. ( f44| ) that 



exp(a — j)x = cosh \[2x H — sinh v^x. (358) 
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Substituting in the relation exp(ce— 7)2; = exp ax exp(— 72;) the expression of the exponentials 
from Eqs. (|337|), (|338|) and fl358|) yields 



/In + /II + /I2 + /Is = ^sh ^2; 



/40/41 - /40/43 + /41/42 + /42/43 = -y= sinh \/2x, 



(359) 
(360) 



where /40, /41, /42, /43 are functions of x. From relations (|359j ) and ( p60|) it can be inferred 
that 



f . hi - /43 \ 2 . / f . hi + /43 \ 2 

/40 H 7= — + /42 H t=f — = exp V2x, 



V2 



/4i-/- 2 



y/2 



10 



43 \ . / , /41 + /43 \ / nz n 

7H + v /42 ~^/H = ex P("^)' 



which means that 

'/40 + 



/41 — /43 X 



+ /42 + 



/.I • / 



/ 13 



40 



/41 — /43 X 



V2 



+ / 



'42 



(361) 
(362) 



/41 + /43 N 



V2 



1.(363) 



An equivalent form of Eq. ( 363 ) is 

/40+/41+/42+/43+ 2 (/Io/l2+/Il/l3)+ 4 (/Io/4l/43+/4o/42/l3-/4o/Il/42-/4l/l2/43) = 1.(364) 



The form of this relation is similar to the expression in Eq. (303). Similarly, since 

(a + 7 ) 2m = (-l) m 2 m , (a + 7 ) 2m+1 = (-l) m 2 m (a + 7), (365) 



(366) 



it can be shown from the definition of the exponential function, Eq. ( |44"| ) that 
exp(a + 7)2; = cos V2x H — sin V2x. 



Substituting in the relation exp(a + 7)2; = exp ax exp 72; the expression of the exponentials 
from Eqs. ( |337| ) , ( |338| ) and ( |366|) yields 

(367) 
(368) 



/lo ~~ /Ii + /I2 ~~ /I3 ~~ cos V%x, 



/40/41 + /40/43 - /41/42 + /42/43 = -/= sin y/2x, 



where /40, /41, fw, /43 are functions of x. 
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Expressions of the four-dimensional cosexponential functions ( p39| )-( p42| ) can be obtained 
using the fact that [(1 + i)/V2] 4 = -1, so that 

fAo( x ) = _ (cosh — ~~r=^x + cos — -r=-%) , (369) 



2 V V2 y/2 

fil (x) = (sinh ^x + sin l -^ X ) , (370) 

/42 (,) = Ifcoshi±!,-cos^V (371) 



2i V ' y/2 \/2 
/ 43 (x) = (sinh ±+± x - sin i±^x ) . (372) 



Using the addition theorems for the functions in the right-hand sides of Eqs. ( |369| )-( |37"2"D , 
the expressions of the four-dimensional cosexponential functions become 

/ 4 o (x) = cos cosh , (373) 

fii(x) = —= | sin —= cosh —= + sinh —= cos — = ) , (374) 

f 42 [x) = sin sinh -^=, (375) 

1 / . X X X X \ 
I sin msn sinn rns 



U 3 (x) = -j= (^sin -j= cosh -j= - sinh -j= cos -j= ) . (376) 

It is remarkable that the series in Eqs. (373)- ( [3761) , in which the terms are either of the 
form x 4m , or x 4m+1 , or x 4m+2 , or x 4m+3 can be expressed in terms of elementary functions 
whose power series are not subject to such restrictions. The graphs of the four-dimensional 
cosexponential functions are shown in Fig. 3. 

It can be checked that the cosexponential functions are solutions of the fourth-order 
differential equation 

d 4 C 

whose solutions are of the form ((u) = Af^(u) + Bfn(u) + Cf 42(11) + Df 43 (u). It can also 
be checked that the derivatives of the cosexponential functions are related by 

d/40 _ f dfa _ , _ , #43 _ , , 970 x 

~i — _ J"43> —, — — 740) — — 741) —, — — 742- 
du du du du 
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4.4 The exponential and trigonometric forms of planar four- 
complex numbers 

Any planar fourcomplex number u = x + ay + (3z + jt can be writen in the form 



x + ay + (3z + jt = e 



x 1 +ayi+f3z 1 +~{t 1 



(379) 



The expressions of xi, y%, Z\,t\ as functions of x, y, z, t can be obtained by developing e ayi , e@ Zl 
and e 7 * 1 with the aid of Eqs. (336)-( p38| ), by multiplying these expressions and separating 
the hypercomplex components, and then substituting the expressions of the four-dimensional 
cosexponential functions, Eqs. fl373|)-(|376|), 



Xl i yi+h yi - h . . yi+h yi-ti 
x = e 1 cos z\ cos = — cosn = sm z\ sin j= — sinn 



V2 y/2 



V2 V2 



(380) 



y = e 



. , , Vi + h . ,yi-h ( 7r\ . yi + h , 
sm \ z\ + — cos = — smn = V cos \ z\ + — sm -= — cosn = — 



Xl i . yi + h yi-h yi + h , yi - h 
z = e 1 cos Z\ sm — — smn -= h sm z\ cos 7= — cosn 



y/2 y/2 



V2 y/2 



(382) 



t = e 



j-i 



, tt\ 2/1 + *i . , yi - h . ( 7r\ . yi + ti yi-h 

cos \ z\ + — cos = — smn = h sm \ z\ + — sm = — cosn = — 

1 ' * * V 1./ y2 s/2 



47 y/2 y/2 
The relations (|380|)-(|3"83|) can be rewritten as 



x + 



z + 



X 



y-t 
y + t 
y-t 
y + t 

V2 



e Xl cos I Z\ + 



e Xl sin I z\ + 



V2 



e *i s - m l Zl -yi±±l) e ~(yi-ti)/V2 



V2 



(384) 
(385) 
(386) 
(387) 



By multiplying the sum of the squares of the first two and of the last two relations ( 384Q -(^87 
it results that 



(388) 



or 



e Xl = p. 



(389) 
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By summing the squares of all relations ( p84| )-( p87| ) it results that 
d 2 = (? cosh \/2(yi —t\) ■ 



(390) 



Then the quantities yi,z±,ti can be expressed in terms of the angles (f>, x,ij) defined in Eqs. 
d3T0D-(p3D as 



V2 



. yi + h 

e (j/i-ti)/V2 
\/2 [cosh y/2{y\ - t x ) 



1/2 



sin?/;, 



-(j/l-tO/V2 



V2 [coshV2(yi -h) 



1/2 



cos V 1 - 



(391) 
(392) 

(393) 



From Eq. (|393| ) it results that 
1 



yi — t\ = —= In tan ip, 
v2 



so that 
Vl 



+ x . i , , , <^-x , <A + X 
+ - — m tan ^, z\ = — - — , t\ 



In tan ^. 



2^2 2^/2 2 ' 2^/2 2^2' 

Substituting the expressions of the quantities x\, y\, Zx,t% in Eq. (|379[) yields 



u = p exp 



1 1 / q + 7 

27=(a- 7 )lntan^ + -(^ + - 7 =- 



1 fa a + 7 , 



(394) 



(395) 



(396) 



which will be called the exponential form of the planar fourcomplex number u. It can be 
checked that 

1 a — 7 fl a — 7 \ / /3 a + 7\ 



exp 



exp 



1(b + 2±1 



1 fa a + 7 , 



2^2 + U + cos * H 2 + ivr ) si " *■ (397) 



1 a — 7 /l a — 7 \ 



: + - 



■ + 



2 ' 2^2 ' V2 2^/2 / V2 2^2 

which shows that e [P+(^+j)/V2\<l>/2 and e -[/3-(a+ 7 )/ v / 2]x/2 

are periodic functions of (f) and 

respectively x, with period 2tt. 



cosx- 



/3 a + 7 



sin x, (398) 



The exponential of the logarithmic term in Eq. ( 396 ) can be expanded with the aid of 
the relation (|35Sj ) as 



exp 



2^2 



(a — 7) In tan ip 



(sin 2^) 1/2 



, 7T \ a — 7 . / 7T 

cos(^--j+^ r sm^-^ 



. (399) 
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Since according to Eq. ( |32C| ) tan = p + /p_, then 

. , , P+P- 
sin w cos w = k- , 

p 2 , + p 2 ' 
and it can be checked that 



(400) 



p 2 + + pl =2d\ 



(401) 



where d has been defined in Eq. ( |308| ) . Thus 



p 2 = d 2 sin 20, 



(402) 



so that the planar fourcomplex number u can be written as 



u = d 



, n \ a - 7 ■ / , 7r 
co8 (^_-j+ sin (^-- 



exp 



1 a + 7\, 1 a + 7 



(403) 



which will be called the trigonometric form of the planar fourcomplex number u. 

If ui,U2 are planar fourcomplex numbers of moduli and angles d\, cfti, xi, 01 and re- 
spectively d2,4>2,X2,ip2, the product of the factors depending on the planar angles can be 
calculated to be 



OL — ^ OL — ^ 

[cos(0i - 7T/4) H — -j=- sin(0i - 7r/4)][cos(0 2 - tt/4) + -j= sin(0 2 - ^/4)] 



[cos(0l - 2 



a — 7 



COs(0l + 2 )]. 



The right-hand side of Eq. (f404|) can be written as 



(404) 



ol — *y 

COs(0l - 02) 7=- COS (-01 + 02) 

A/2 

= [2(cos 2 0i cos 2 02 + sin 2 0i sin 2 2 )] 1/2 [cos(0 - vr/4) + - 7 sin(0 - 7r/4)](405) 



where the angle 0, determined by the condition that 



tan(0 — 7r/4) = — cos(-0i + "02)/ cos(0i — 02 ) 



(406) 



is given by tan0 = tan0i tan 02 , which is consistent with Eq. ( ^22|) . The modulus d of the 
product u\U2 is then 



d = \Pld\d2 (cos 2 0i cos 2 02 + sin 2 0i sin 2 02 



\l/2 



(407) 
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4.5 Elementary functions of planar fourcomplex variables 

The logarithm u% of the planar fourcomplex number u, U\ = Inn, can be defined as the 
solution of the equation 

u = e Ul , (408) 



written explicitly previously in Eq. (|37S| ), for u\ as a function of u. From Eq. ( |396| ) it results 
that 

\nu = \n P +^(a- 1 )\nt^+-{^ + —j^^--^3-—j^Jx, (409) 
which is multivalued because of the presence of the terms proportional to <ft and x- It can 



be inferred from Eqs. (321) and ( 322 ) that 



ln(im') = Inu + Inn', (410) 

up to multiples of n[/3 + (a + j)/V2\ and n[fl - (a + j)/V^]- 
The power function u m can be defined for real values of n as 

u m = e mlnu . (411) 

The power function is multivalued unless n is an integer. For integer n, it can be inferred 
from Eq. (ElC]) that 



(uu') m = u m u' m . (412) 

If, for example, m = 2, it can be checked with the aid of Eq. ( |403[) that Eq. ( |411[) gives 
indeed (x + ay + j3z + ~ft) 2 = x 2 - z 2 - 2yt + 2a(xy - zt) + (5{y 2 - t 2 + 2xz) + 27(2^ + yz). 

The trigonometric functions of the fourcomplex variable u and the addition theorems for 
these functions have been written in Eqs. (|90|)-(p^). The cosine and sine functions of the 
hypercomplex variables ay, (3z and jt can be expressed as 

cos ay = f 40 (y) - (3f A2 {y), sin ay = a/ 4 i(y) - 7/43(2/), (413) 
cosfiz = coshz, s\n[3z = /?sinhz, (414) 
cos jt = f m (t) + j3 f 42 (t) , sin 7 t = 7/41 (t) - a/43 (*) ■ (415) 
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The cosine and sine functions of a planar fourcomplex number x + ay + f3z + jt can then be 
expressed in terms of elementary functions with the aid of the addition theorems Eqs. (92), 
(|9"3|) and of the expressions in Eqs. (f413|)-(f4l"5|). 



The hyperbolic functions of the fourcomplex variable u and the addition theorems for 



these functions have been written in Eqs. (M)-(IOO). The hyperbolic cosine and sine func- 



tions of the hypercomplex variables ay, f3z and yt can be expressed as 

cosh ay = / 40 (y ) + /3f 42 (y) , sinh ay = af 41 (y) + 7/43 {y ) , (416) 

cosh/3z = cosz, sinh/3z = (3s'mz, (417) 

cosh jt = / 40 (i) - /3/ 42 (t) , sinh jt = 7/41 (t) + a/43 (t ) • (418) 

The hyperbolic cosine and sine functions of a planar fourcomplex number x + ay + f3z + jt 
can then be expressed in terms of elementary functions with the aid of the addition theorems 
Eqs. (H), ( jjjXj ) and of the expressions in Eqs. (gl^)-(|5g). 

4.6 Power series of planar fourcomplex variables 

A planar fourcomplex series is an infinite sum of the form 

a + a\ + a 2 H h aH , (419) 



where the coefficients a; are planar fourcomplex numbers. The convergence of the series (|419|) 
can be defined in terms of the convergence of its 4 real components. The convergence of a 
planar fourcomplex series can however be studied using planar fourcomplex variables. The 
main criterion for absolute convergence remains the comparison theorem, but this requires a 
number of inequalities which will be discussed further. 

The modulus of a planar fourcomplex number u = x + ay + (3z + jt can be defined as 

| M | = (x 2 +y 2 + z2+t2) l/ 2) (42Q) 

so that, according to Eq. ( |308j ) , d = \u\. Since \x\ < < < l n M^I — M> a 

property of absolute convergence established via a comparison theorem based on the modulus 



of the series ( 419 ) will ensure the absolute convergence of each real component of that series. 
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The modulus of the sum u± + u 2 of the planar fourcomplex numbers u\,u 2 fulfils the 
inequality 

\\u\\ - 1 1 1 < \u\ + u 2 \ < \ui\ + \u 2 \. (421) 
For the product the relation is 

\u lU2 \ < V2\ Ul \\u 2 \, (422) 



as can be shown from Eq. ( |40 7j ) . The relation ( (422|) replaces the relation of equality extant 
for regular complex numbers. The equality in Eq. (fl22j ) takes place for cos 2 (V'i — ip 2 ) = 
1, cos 2 ('0i + ip 2 ) = 1) which means that x\ + (yi — ti)/\^2 = 0, Z\ + (yi + t\)/^/2 = 0, x 2 + 
(2/2 - h)/V2 = 0, z 2 + (y 2 + t 2 )/V2 = 0, or Xl - {y x - h)/V2 = 0, Zl - ( Vl + h)/V2 = 0,x 2 - 
(lJ2 — ^2)/ V% = 0, z 2 — (y 2 + t 2 )j \f2 = 0. The modulus of the product, which has the property 
that < |«iit2|) becomes equal to zero for cos 2 ^ — ip 2 ) = 0; cos 2 ^ +^2) = 3 which means 
that xi + { yi -h)/V2 = 0,zt + (yi+h)/V2 = 0,x 2 - {y 2 - 1 2 )/^2 = 0,z 2 - (y 2 + t 2 )/V2 = 0, 
or x x - ( yi - h)/V2 = 0, Zl - (yi + t{)/y/2 = 0,x 2 + (y 2 - t 2 )/V2 = 0,z 2 + (y 2 + t 2 )/y/2 = 0. 
as discussed after Eq. (|307| ). 
It can be shown that 

x 2 + y 2 + z 2 + t 2 < \u 2 \ < V2(x 2 + y 2 + z 2 + t 2 ). (423) 

The left relation in Eq. (|423| ) becomes an equality for sin 2 2^ = 1, when p + = which 
means that x(y — t) + z(y + 1) = 0. The right relation in Eq. (142$) becomes an equality for 
sin 2 2-0 = 0, when x+(y-t)/y/2 = 0, z+(y+t)/y/2 = 0, or x-(y-t)/y/2 = 0, z-(y+t)/y/2 = 
0. The inequality in Eq. (|422j ) implies that 

|u'| < 2 ( ^ 1)/2 |u| / . (424) 



From Eqs. (|422[) and fl424| ) it results that 

\au l \ <2 l/2 \a\\u\ l . (425) 
A power series of the planar fourcomplex variable u is a series of the form 

ao + a\u + a 2 u 2 + • • • + aiv! + • • • . (426) 
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Since 



1=0 



< 



E 2V2 i 



ai\\u\ 



1=0 



a sufficient condition for the absolute convergence of this series is that 

\/2|ai+i||«| 1 
hm j — j < 1. 

Thus the series is absolutely convergent for 



(427) 



(428) 



\u\ < c, 



(429) 



where 



lim 



\ a l\ 



(430) 



The convergence of the series ( f426| ) can be also studied with the aid of the transformation 



x 



+ ay + f3z + jt = V2(ei£ + e\v + e 2 r + e 2 C), 



(431) 



where £, v, r, £ have been defined in Eq. (|309|) , and 



1 a — 7 
ei = 2 + W 61 

It can be checked that 



P a + 7 1 a — 7 

2 + ^7T 62 = 2 ~ ^7T' 



C2 



^ _l_ a + 7 
'2 2^2 ' 



(432) 



2 ~2 

e x = ex, e 1 



-ei, eiei = ei, e| = e 2 , e| 



-e 2 , e 2 e 2 = e 2 , 



eie 2 = 0, eie 2 = 0, eie 2 = 0, e 2 ei = 0. 



(433) 



The moduli of the bases in Eq. ( |432| ) are 



1 



1 



1 



e i 



y/2' 1611 V2' 
and it can be checked that 



e 2 



e 2 



1 

71' 



(434) 



\x + ay + pz + jt\ 2 =f + v 2 + T 2 + C 2 . 



(435) 



The ensemble ei,ei,e 2 ,e 2 will be called the canonical planar fourcomplex base, and Eq. 



(431) gives the canonical form of the planar fourcomplex number. 

If u = u'u" ', the components £, v,t, £ are related, according to Eqs. ( 325[ )-( [32"8| ) by 



£ = ^-„V),« = x/2(£V' WO, r = ^(rV'-C'C"), C = \/2(r'C"+Cr), (436) 
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which show that, upon multiplication, the components £, v and r, £ obey, up to a normaliza- 
tion constant, the same rules as the real and imaginary components of usual, two-dimensional 
complex numbers. 



If the coefficients in Eq. (426) are 
ai = m + aa tl + fia i2 + 7^3, 



and 



„ ill - a l3 ~ an + a i3 a n - a l3 ~ 
An = ai -\ 7= — , An = ai2-\ 7= — , A i2 = a i0 7= — , A i2 



V2 

the series (426) can be written as 



V2 



V2 



(437) 



. a ll + a l3 / AQQ\ 

-a«2H 7= — , (438) 



2 i/2 [(eiAn + eil a )(ei£ + eiv) 1 + {e 2 A l2 + e 2 A l2 )(e 2 T + e 2 Q l 

1=0 



Thus, the series in Eqs. ([426; ) and ( [439|) are absolutely convergent for 

P+ < ci, p- < c 2 , 

where 



ci = lim 



A 2 + A 2 



1/2 



A 2 +l l + Af +l l 



A 2 + A 2 

A 12 * 



12 



1/2 



1/2' 



V2 



(439) 



(440) 



(441) 



^4 2 +i,2 + Af + i j2 

It can be shown that c = (l/v2)min(ci, C2), where min designates the smallest of the 
numbers ci,c 2 . Using the expression of |tt| in Eq. ( |435| ), it can be seen that the spherical 
region of convergence defined in Eqs. ( |429j ), ( |430| ) is included in the cylindrical region of 
convergence defined in Eqs. ( |44C| ) and ( [44 1| ). 



4.7 Analytic functions of planar fourcomplex variables 

The fourcomplex function f(u) of the fourcomplex variable u has been expressed in Eq. ( |136j ) 

in terms of the real functions P(x, y, z, t), Q(x, y, z, t), R(x, y, z, t), S(x, y, z, t) of real variables 

x, y, z, t. The relations between the partial derivatives of the functions P, Q, R, S are obtained 

by setting succesively in Eq. (|137[ ) Ax — > 0, Ay = Az = At = 0; then Ay — > 0, Ax = Az = 

At = 0; then Az -> 0, Ax = Ay = At = 0; and finally At -> 0, Ax = Ay = Az = 0. The 

relations are 

dP = dQ = dR = dS 
dx dy dz dt ' 

63 



The relations (|442|)-(p5D are analogous to the Riemann relations for the real and imagi- 



nary components of a complex function. It can be shown from Eqs. (442)-(445) that the 
component P is a solution of the equations 

Q2p Q2p Q2p Q2p 

^ + 1^ =0 < w + lv =0 > (446) 

and the components Q, R, S are solutions of similar equations. As can be seen from Eqs. 

(446), the components P, Q, R, S of an analytic function of planar fourcomplex variable are 

harmonic with respect to the pairs of variables x, y and z, t. The component P is also a 

solution of the mixed-derivative equations 

d 2 P _ _(PP_ d^P _ d 2 P d 2 P _ c?P_ (PP_ _ d 2 P 
dx 2 dydt 1 dy 2 dxdz' dz 2 dydt' dt 2 dxdz' 

and the components Q, R, S are solutions of similar equations. The component P is also a 

solution of the mixed-derivative equations 

d 2 P d 2 P d 2 P _ d 2 P 
dxdy dzdV dxdt dydz' 

and the components Q, R, S are solutions of similar equations. 



4.8 Integrals of functions of planar fourcomplex variables 

The singularities of planar fourcomplex functions arise from terms of the form l/(u — uo) m , 
with m > 0. Functions containing such terms are singular not only at u = Uq, but also at all 
points of the two-dimensional hyperplanes passing through uq and which are parallel to the 
nodal hyperplanes. 

The integral of a planar fourcomplex function between two points A, B along a path 
situated in a region free of singularities is independent of path, which means that the integral 
of an analytic function along a loop situated in a region free from singularities is zero, 

j f{u)du = 0, (449) 

64 



where it is supposed that a surface E spanning the closed loop T is not intersected by any 
of the two-dimensional hyperplanes associated with the singularities of the function f(u). 
Using the expression, Eq. ( pL36| ) for f(u) and the fact that du = dx + ady + f3dz + 'ydt, the 
explicit form of the integral in Eq. ( |44S| ) is 

j> f(u)du = j> [(Pdx - Sdy - Rdz - Qdt) + a(Qdx + Pdy - Sdz - Rdt) 

+(3(Rdx + Qdy + Pdz - Sdt) + j(Sdx + Rdy + Qdz + Pdt)} . (450) 

If the functions P,Q,R,S are regular on a surface X spanning the loop I, the integral along 
the loop r can be transformed with the aid of the theorem of Stokes in an integral over the 
surface £ of terms of the form dP/dy + dS/dx, dP/dz + dR/dx, dP/dt + dQ/dx, dR/dy — 
dS/dz, dS/dt — dQ/dy, dR/dt — dQ/dz and of similar terms arising from the a,f3 and 7 
components, which are equal to zero by Eqs. (|442|)-(^45|), and this proves Eq. ([449|). 



The integral of the function (u — uq)" 1 on a closed loop T is equal to zero for m a positive 
or negative integer not equal to -1, 

<j> (u — uo) m du = 0, m integer, m 7^ —1. (451) 

This is due to the fact that J(u — uo) m du = (u — uo) m+l /(m + 1), and to the fact that the 
function (u — uo) m+1 is singlevalued for m an integer. 

The integral § du/{u — uq) can be calculated using the exponential form ( |396| ), 



u — uq = p exp 
so that 



du dp 1 l/o a + 7\ 1 / a + 7 s 



(452) 



u- „„ p + 2vi (a " 7)dlntm * + 5 (" + ^r) * - 5 (" - ^r) (453> 

Since p and V 1 are singlevalued variables, it follows that j> r dp/ p = 0, j> r din tan ijj = 0. On the 
other hand, 4> an d X are cyclic variables, so that they may give a contribution to the integral 
around the closed loop V. Thus, if C+ is a circle of radius r parallel to the £Ov plane, whose 
projection of the center of this circle on the t;Ov plane coincides with the projection of the 
point no on this plane, the points of the circle C + are described according to Eqs. (|309|) -( ^13|) 
by the equations 

£ = £0 + r sin ip cos </>, v = v + r sin ip sin </>, t = tq + r cos ip cos Xi 

C = Co + rcosipsmx, (454) 
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for constant values of % an d ip, ifi 7^ 0, 7r/2, where uo = xo + ayo + (3zq + 7*0, and £o> v 0, to, Co 
are calculated from xo,yo 5 -20)*o according to Eqs. (|30S|) . Then 



c+u-uq 



7T /? + 



a + 7 



(455) 



If C_ is a circle of radius r parallel to the rO£ plane, whose projection of the center of 
this circle on the tOC, plane coincides with the projection of the point uq on this plane, the 
points of the circle C_ are described by the same Eqs. (|454| ) but for constant values of (j> 
and tp, ijj 7^ 0,7r/2. Then 

(in / „ a + 7 " 



-7T /? 



(456) 



/c_ u — uo \ \f2 

The expression of fpdu/(u — uq) can be written as a single equation with the aid of the 



functional int(M, C) defined in Eq. (153) as 



du 



r u - u 



vr [(3 + 



a + 7 



int(uo^,r 5 „) -ir (ft 



a + 7 



int(uorC: r Tc)> (457) 



where Uoft>, UotC and T^, are respectively the projections of the point uq and of the loop 
r on the planes and r£. 

If /(u) is an analytic planar fourcomplex function which can be expanded in a series as 



written in Eq. (131), and the expansion holds on the curve T and on a surface spanning T, 
then from Eqs. (|45l| ) and (|457|) it follows that 



f(u)du 

T U-Uq 



7T 



a + 7 



int(ito£«,IV 







a + 7 



int(norC 5 r rc) 



/M,(458) 



where r^„,r r( ^ are the projections of the curve V on the planes and respectively t£, as 
shown in Fig. 2. As remarked previously, the definition of the variables in this section is 
different from the former definition for the circular hypercomplex numbers. 

Substituting in the right-hand side of Eq. (|45§| ) the expression of f(u) in terms of the 



real components P,Q,R,S, Eq. (|136[) , yields 
f{u)du 



r u - Uq 



7T 



-7T 



P + 



a + 7 
a + 7' 



P 
P- 



1 + 



a — 7 
a — 7 



R 
R- 



7 

7" 



1-/9 

1 + P 
V2 



Q 
Q- 



a 



a 



1 + /3 

l-P 
V2 



int(w ^,r 

int(n 0T c,r r 
(4 
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where P, Q, R, S are the values of the components of / at it = uq. 



If f(u) can be expanded as written in Eq. (131) on T and on a surface spanning T, then 
from Eqs. ( |451| ) and (457) it also results that 



f(u)du TT 



{u — uo) m+1 ml 



(3 H -=- int(u ^, T^) -1/3 ) int(-u 0r c, T tC ) 



V2 J V" ^2 

where it has been used the fact that the derivative f^ m \uo) of order n of f(u) at u = uq is 
related to the expansion coefficient in Eq. (|l3lj ) according to Eq. (|135| ). 

If a function f(u) is expanded in positive and negative powers of u — Uj, where Uj are 
planar fourcomplex constants, j being an index, the integral of / on a closed loop T is 
determined by the terms in the expansion of / which are of the form aj/(u — uj), 



/<•">(«„), (460) 



3 J 

Then the integral of / on a closed loop T is 

j>^f{u)du = tt(p+ "J^ 7 ^ ^int(u j?l ,,r^)aj-7r - "J^ 7 ^ ^] int(-u iTC , r T? )aj.(462) 

4.9 Factorization of planar fourcomplex polynomials 

A polynomial of degree m of the planar fourcomplex variable u = x + ay + (3z + yt has the 
form 

P m (u) = u m + a lU m ~ l + ■■■ + a m -iu + a m , (463) 

where the constants are in general planar fourcomplex numbers. 

It can be shown that any planar fourcomplex polynomial has a planar fourcomplex root, 
whence it follows that a polynomial of degree m can be written as a product of m linear 
factors of the form u — uj, where the planar fourcomplex numbers Uj are the roots of the 
polynomials, although the factorization may not be unique, 

in 

Pm(u) = Y[(u- Uj ). (464) 
3=1 

The fact that any planar fourcomplex polynomial has a root can be shown by considering 
the transformation of a fourdimensional sphere with the center at the origin by the function 
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u m . The points of the hypersphere of radius d are of the form written in Eq. (403), with d 
constant and (f>,x,ift arbitrary. The point u m is 

q — 7 



u m = d ra 



7T 



cos(^--j+ ^ 



• sin if) 



7T 



exp 



1 ( n " + 7\ , 1 ( n « + 7 



It can be shown with the aid of Eq. ( 407 ) that 
uexp 



2V + V2 



1 fa a + ^ , 



(465) 



(466) 



so that 



cos(V> - 7r/4) + - — -(if) - it/A) 
sin 



cxp 



m(f)--[P 



a + 7 
~72 



a — 7 

cos(if) — vr/4) H sin(^ — tt/4) 

v2 



The right-hand side of Eq. ( [46 7|) is 

2 » 



a — 7 

cos e H =- sin e 

V2 



^C^cos" 



■ 2fc 

esm e, 



(467) 



(468) 



fc=0 
i2k \ /^ffc 



where e = ^ — 7r/4, and since C^m — C„, it can be concluded that 

2 



a — 7 
cos e H t=- sm e 



> 1. 



(469) 



Then 



d m < \u m \ < 2^ m - 1 ^ 2 d m , 



(470) 



which shows that the image of a four-dimensional sphere via the transformation operated by 
the function u m is a finite hypersurface. 
If u' = u m , and 



u = d 



(x — 7 

cos(if)' — vr/4) H t=— sin('0 / — 7r/4) 

v2 



exp 



2T 



then 



(471) 



cj)' = m<p, x = m Xi t an V* = t an?Tl V'- 



(472) 



Since for any values of the angles (fi',x',ip' there is a set of solutions (f), x, V 1 °f Eqs. ( f472| ), 
and since the image of the hypersphere is a finite hypersurface, it follows that the image of 
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the four-dimensional sphere via the function u m is also a closed hypersurface. A continuous 
hypersurface is called closed when any ray issued from the origin intersects that surface at 
least once in the finite part of the space. 

A transformation of the four-dimensional space by the polynomial P m {u) will be con- 
sidered further. By this transformation, a hypersphere of radius d having the center at the 
origin is changed into a certain finite closed surface, as discussed previously. The transfor- 
mation of the four-dimensional space by the polynomial P m (u) associates to the point u = 
the point /(0) = a m , and the image of a hypersphere of very large radius d can be repre- 
sented with good approximation by the image of that hypersphere by the function u m . The 
origin of the axes is an inner point of the latter image. If the radius of the hypersphere is 
now reduced continuously from the initial very large values to zero, the image hypersphere 
encloses initially the origin, but the image shrinks to a m when the radius approaches the 
value zero. Thus, the origin is initially inside the image hypersurface, and it lies outside the 
image hypersurface when the radius of the hypersphere tends to zero. Then since the image 
hypersurface is closed, the image surface must intersect at some stage the origin of the axes, 



which means that there is a point u\ such that f(u\) = 0. The factorization in Eq. (464) 
can then be obtained by iterations. 

The roots of the polynomial P m can be obtained by the following method. If the constants 
in Eq. ( j463| ) are a/ = a/o + cra/i + f3ai 2 + 70/3, and with the notations of Eq. fl438p , the 
polynomial P m (u) can be written as 



\m—l 



P m = Yl 2 (m -' )/2 ( ei Ai + giin)(ei£ + 
1=0 

m 

+ ^ 2 (m-0/2 (e2 ^ 2 + e 2 A l2 )(e 2 T + e 2 () m ~ l , (473) 
1=0 

where the constants A[k,Aik,k = 1,2 are real numbers. Each of the polynomials of degree 
m in ei£ + e\v, e 2 T + e 2 Q in Eq. fl473j ) can always be written as a product of linear factors of 
the form ei(£ — £ p ) + e\(v — v p ) and respectively e 2 {r — T p ) + e 2 (C ~ Cp)> where the constants 
Cp,v P ,Tp,Cp are real, 

m m 

£ 2^l\e l A ll + eiiu)( ei £ + ^v) m ~ l = J] 2 m/2 {e x (£ - + e 1 (v - v p )} , (474) 
1=0 p=i 
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J2 2^/ 2 (e 2 A l2 + e 2 A l2 )(e 2 r + ~e 2 Q m ~ l = J[ 2 m ' 2 {e 2 (r - r p ) + e 2 (( - ( p )} . (475) 

1=0 p=l 

Due to the relations ( |433| ), the polynomial P m {u) can be written as a product of factors 
of the form 



P m (u) = J] 2 m/2 {ei(£ - Z P ) +h(v- v P ) + e 2 (r - r p ) + e 2 (C - C P )} • 
p=i 

This relation can be written with the aid of Eq. (431) in the form (464), where 



(476) 



'2{eii p + hvp + e 2 T p + e 2 C P )- (477) 

The roots ei^ p + eif p and e 2 r v + e 2 C p defined in Eqs. ( |474| ) and respectively (|475| ) may 
be ordered arbitrarily. This means that Eq. ( (477|) gives sets of m roots u\,...,u m of the 
polynomial P m (tt), corresponding to the various ways in which the roots ei£ p + e\v p and 
&2 T p + e 2 Cp are ordered according to p for each polynomial. Thus, while the hypercomplex 
components in Eqs. ( |474| ) , Q475 ) taken separately have unique factorizations, the polynomial 
P m (u) can be written in many different ways as a product of linear factors. The result of 
the planar fourcomplex integration, Eq. Q462 ), is however unique. 

If, for example, P(u) = u 2 + l, the possible factorizations are P = (u— e\ — e 2 )(u + ei + e 2 ) 
and P = (u — e\ + e 2 )(u + e\ — e 2 ) which can also be written as u 2 + 1 = (u — (3)(u + (3) 
or as v? + 1 = ju — (a + 7)/v / 2j |(u + (a + r y)/^/2^. The result of the planar fourcomplex 



integration, Eq. ( [462|) , is however unique. It can be checked that (±ei±e2 



7, \2 



-e\— e 2 



4.10 Representation of planar fourcomplex numbers by irre- 
ducible matrices 



If T is the unitary matrix, 





( i 


1 

2 





T = 




1 


1 

2 

1 

2 


i 






{ 


1 

2 


1 

V2 


can be shown that the matr 


TUT 




V 


■ 

^2 / 



1 \ 

2 



(478) 



(479) 
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where U is the matrix in Eq. ( |332| ) used to represent the planar fourcomplex number it. In 
Eq. ( 479 ) , Vi , V2 are the matrices 



Vi 



f 1 v-t 
X+ V2 



z + 



y+t 
V2 



V, 



V 



-z 



V2 X ^ V2 



v-t 

y+t 
V2 



r — £=* 
E V2 



(480) 



In Eq. ( |479| ) , the symbols denote the matrix 




(481) 



The relations between the variables x+(y—t)/y/2, z+{y+t)/y/2, x—(y—t)/\ / 2, —z+(y+t)/\ / 2 
for the multiplication of planar fourcomplex numbers have been written in Eqs. ( |325| )-( j32"8|) . 
The matrix TUT -1 provides an irreducible representation Q of the planar fourcomplex 
number u in terms of matrices with real coefficients. 



5 Polar omplex Numbers in Four Dimensions 
5.1 Operations with polar fourcomplex numbers 

A polar fourcomplex number is determined by its four components (x,y,z,t). The sum of 
the polar fourcomplex numbers (x,y,z,t) and (x', y' , z 1 , if) is the polar fourcomplex number 
(x + x' , y + y' , z + z' ,t + t'). The product of the polar fourcomplex numbers (x, y, z, t) and 
(a/, y 1 , z', t') is defined in this work to be the polar fourcomplex number (xx' + yi! + zz 1 + 
ty' , xy' + yx' + zt' + tz' , xz' + yy' + zx' + tt', xt' + yz' + zy' + tx'). Polar fourcomplex numbers 
and their operations can be represented by writing the polar fourcomplex number (x, y, z, t) 
as u = x + ay + (3z + jt, where a, (3 and 7 are bases for which the multiplication rules are 

a 2 = p, (3 2 = 1, 7 2 = (3,a(3 = (3a = 7, 07 = 7a = -1, /3 7 = 7 /3 = a. (482) 

Two polar fourcomplex numbers u = x+ay+flz+jt, u' = x' +ay' +(5z' +^t' are equal, u = u', 
if and only if x = x', y = y', z = z' ,t = t'. If u = x + ay + (3z + jt,u' = x' + ay' + f3 'z' + jt' are 
polar fourcomplex numbers, the sum u+u' and the product uvl defined above can be obtained 
by applying the usual algebraic rules to the sum (x + ay + j3z + jt) + (x' + ay' + (3z' + jt') 
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and to the product [x + ay + (3z + "ft)(x' + ay' + fiz' + "ft'), and grouping of the resulting 
terms, 

u + u ' = x + x' + a(y + y') + 0{z + z) + -f(t + t'), (483) 



uu = xx + yt' + zz' + ty' + a{xy + yx + zt' + te') + f3(xz + yy' + zx + tt!) 

+"f(xt' + yz + zy' + tx). (484) 

If u,u',u" are polar fourcomplex numbers, the multiplication is associative 

(uu')u" = u{u'u") (485) 
and commutative 

uu = uu, (486) 

as can be checked through direct calculation. The polar fourcomplex zero is 0+a-0+/3-0+7-0, 
denoted simply 0, and the polar fourcomplex unity is l + a- + /3-0 + 7-0, denoted simply 
1. 

The inverse of the polar fourcomplex number u = x + ay + f3z + "ft is a polar fourcomplex 
number u' = x' + ay' + (3z' + "ft' having the property that 

uu' = 1. (487) 

Written on components, the condition, Eq. ([48 7|), is 



(488) 



xx' + ty' + zz' + yt' = 1, 
yx' + xy' + tz' + zt' = 0, 
zx' + yy' + xz' + tt' = 0, 
te' + zy' + yz' + xt' = 0. 

The system ( |488f ) has the solution 

x , = ^(x 2 - z 2 ) + z(y 2 + t 2 ) - 2xyt (489) 

y , = -y(x 2 + z 2 )+t(y 2 -t 2 ) + 2xzt 
v 

z , = -z(x 2 -z 2 )+x(y 2 + t 2 )-2yzt 



v 
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, _ -t(x 2 + z 2 ) - y(y 2 - r) + 2xyz 



(492) 



provided that v ^ 0, where 

v = x 4 + z 4 - y 4 - t A - 2x 2 z 2 + 2y 2 t 2 - Ax 2 yt - Ayz 2 t + 4xy 2 z + 4xzt 2 . (493) 
The quantity v can be written as 

v = v + v-fjL 2 + , (494) 

where 

v + = x + y + z + t, V- = x — y + z — t, (495) 

and 

VL+ = (x- zf + (y - tf. (496) 

Then a polar fourcomplex number q = x + ay + /3z + 7t has an inverse, unless 

u + = 0, or v- = 0, or fj, + = 0. (497) 

The condition u + = represents the 3-dimensional hyperplane x + y + z + t = 0, the condition 
i>_ = represents the 3-dimensional hyperplane x — y + z — t = 0, and the condition = 
represents the 2-dimensional hyperplane x = z,y = t. For arbitrary values of the variables 
x, y, z, t, the quantity v can be positive or negative. If v > 0, the quantity p = v 1 ^ will be 
called amplitude of the polar fourcomplex number x + ay + (3z + jt. Because of conditions 



(|497|) , these hyperplanes will be called nodal hyperplanes. 

It can be shown that if uv! = then either u = 0, or v! = 0, or the polar fourcomplex 
numbers u, u' belong to different members of the pairs of orthogonal hypersurfaces listed 
further, 

x + y + z + t = and x = y = z = t' , (498) 

x — y + z — t = and x = —y = z = —t. (499) 

Divisors of zero also exist if the polar fourcomplex numbers u, u' belong to different members 
of the pair of two-dimensional hypersurfaces, 

x - z = 0, y - t = and x' + z' = 0, y' + t' = 0. (500) 
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5.2 Geometric representation of polar fourcomplex numbers 

The polar fourcomplex number x + ay + (3z + jt can be represented by the point A of 
coordinates (x, y, z, t). If O is the origin of the four-dimensional space x, y, z, t, the distance 
from A to the origin O can be taken as 

<f = x 2 + y 2 + z 2 + t 2 . (501) 

The distance d will be called modulus of the polar fourcomplex number x + ay + (3z + yt, 
d = \u\. 

liu = x + ay + /3z + jt,ui = xx + ayx +f3z\ + jtx 1 u 2 = x 2 + ay 2 + (3z 2 + jt 2 , and u = uiu 2 , 
and if 

Sj + = xj + yj + zj + tj, sj- = Xj - yj + Zj - tj, (502) 
for j = 1,2, it can be shown that 

V + = S1+S2+, V- = S1-S2-. (503) 

The relations (|503|) are a consequence of the identities 

{x\X2 + ^iz 2 + tiy 2 + 2/1*2) + (xiy 2 + Vix 2 + «i*2 + M2) 

+ (x 1 z 2 + zix 2 + yiy 2 + ht 2 ) + (xit 2 + hx 2 + z x y 2 + yiz 2 ) 

= (xi + yi + zx + ti)(s 2 + 2/2 + z 2 + t 2 ), (504) 

{xiX 2 + Z 1 Z 2 + tiJ/2 + 2/1*2) - OciJ/2 + 2/1^2 + ^1*2 + M2) 

+ (x 1 z 2 + zix 2 + yiy 2 + ht 2 ) - (xit 2 + t\x 2 + z x y 2 + y 1 z 2 ) 

= (xx +zx-yx~ h)(x 2 + z 2 -y 2 - t 2 ). (505) 
The differences 

v\ = x — z, vx = y — t (506) 



can be written with the aid of the radius Eq. ( 496 ), and of the azimuthal angle (f>, where 
< (j) < 2ir as 

vx = [J>+ cos (/), vx = fJ>+ sin 0. (507) 
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The variables v+,V—,vi,vx will be called canonical polar fourcomplex variables. The distance 
d, Eq. fl501| ), can then be written as 

(508) 



It can be shown that if u\ = 27 + ay\ + f3z\ + 7*1,1*2 = x 2 + ay 2 + /?£2 + 7*2 are polar 
fourcomplex numbers of polar radii and angles pi-,(j)\ and respectively p 2 -,(p 2 , then the 
polar radius p and the angle <p of the product polar fourcomplex number u\u 2 are 



p + = P1-P2- 



The relation ( |509| ) is a consequence of the identity 

[(xix 2 + z x z 2 + yi* 2 + hy 2 ) - (xiz 2 + z x x 2 + ym + M2)] 2 

+ [(xiV2 + Vix 2 + zit 2 + t\z 2 ) - (27*2 + hx 2 + z x y 2 + y\z 2 )] 2 
= (x! - z x f + ( Vl - *i) 2 l [(x 2 - z 2 f + (y 2 - t 2 ) 2] 



and the relation ( 510 ) is a consequence of the identities 



(27X2 + Z 1 Z 2 + yi*2 + tlZ/2) - (^1^2 + £iX 2 + 2/12/2 + *1*2) 

= (27 - zi)(x 2 - z 2 ) - (yi - *i)(y 2 - h), 



(509) 



(510) 



(511) 



(512) 



(xtyz + I/1X2 + Zl* 2 + tiZ 2 ) - (Xi* 2 + *1X 2 + Ziy 2 + y\Zl) 

= (yt - *i)(x 2 - z 2 ) + (xi - z l )(y 2 - t 2 ). 



(513) 



A consequence of Eqs. ( 503 ) and ( |509D is that if u = u\u 2 , and Uj = Sjs'-pj- , where 
j = 1, 2, then 



1/ = V\V 2 . 



The angles between the line OA and the v + and respectively V- axes are 



(514) 



(515) 
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where < 9+ < it, < #_ < ir. The variable /u+ can be expressed with the aid of Eq. ( |508| ) 

as 



1 1 



-i 



The coordinates x, y, z, t can then be expressed in terms of the distance d, of the polar angles 
9+ , 9- and of the azimuthal angle (j) as 

u + (tan 8 + + tan ) 1 

X = 9 W; fl t fl + 9 ^+ C ° S ^' (517) 
2V2 tan (/-)_ tan y_ 2 

_ /x + (-tanfl + + tan6L) 1 
^ 2^2 tan 9+ tan 2 M+Sm<? ' 1 ' 

/x+(tan0+ + tan0_) 1 

Z = ^W^ — T~+ — a o/ i + cos ^ ( 519 ) 

2v 2 tan 9 + tan 0_ 2 

u+(— tan 9 + + tan 0_) 1 , 

* = 9 fit a + o ~ 2 ^ + sin ^ (520) 
2y2 tan 0+ tan 0_ 2 



If u = u\U2, then Eqs. ( [50 3| ) and fl509| ) imply that 

1 1 

tan0 + = — ^= tan 0i + tan 02+, tan0_ = — y= tan 9\- tan 02- , (521) 

where 

tan^ + = ^+,tan^_ = ^^. (522) 

Sj Sj 

An alternative choice of the angular variables is 

fi + = Videos 9, v + = 2cisin0cosA, V- = 2dsin0sinA, (523) 
where < 9 < 7r/2, < A < 2ir. If u = U\U2, then 

tan A = tan Ai tanA2, dcos 9 = \f2d\d2 cos 9\ cos 6*2, (524) 

where 

Pj_ = \/2dj cos 0j , Sj = 2dj sin 0j cos Xj , s - = 2dj sin 9j sin Aj , (525) 

for j = 1,2. The coordinates x,y,z,t can then be expressed in terms of the distance d, of 
the polar angles 9X and of the azimuthal angle <\> as 

x = —= sin 9 sin(A + 7r/4) H — ^= cos 6* cos (j), (526) 
V 2 V2 
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y = -^= sin#cos(A + 7r/4) + -^= cos^sin^, (527) 

z = — ^= sin 9 sin(A + vr/4) -= cos cos (f>, (528) 

V2 v2 

t = — sin#cos(A + 7r/4) ^= cos 6* sin (p. (529) 

v 2 v 2 

The polar fourcomplex numbers 

1 + Q + /3 + 7 l-a + /?-7 ,„ n x 
e+ = , e_ = , (530) 

have the property that 

e\ = e+, e?. = e_, e + e_ = 0. (531) 
The polar fourcomplex numbers 

ei = ^— ,ei = ^— (532) 
have the property that 

ef = ei, ef = — ej, eiei = ei. (533) 
The polar fourcomplex numbers e+, e_ are orthogonal to ej, ex, 

e + ex = 0, e + ei = 0, e_ ex = 0, e_ §i = 0. (534) 
The polar fourcomplex number q = x + ay + j3z + jt can then be written as 

x + ay + f3z + jt = v + e + + w_e_ + v\e\ + v\e\. (535) 

The ensemble e+, e_, ei, ei will be called the canonical polar fourcomplex base, and Eq. ( |535| ) 
gives the canonical form of the polar fourcomplex number. Thus, the product of the polar 
fourcomplex numbers u, v! can be expressed as 

uu = v + v' + e + + v-v'_e- + (vxv[ — vxv'x) e i + ( v l^x + v i^l)^ii (536) 

where v ' + = x' + y' + z' + t', v'_ = x' — y' + z' — t' , v[ = x' — y', v[ = z' — t'. The moduli of 
the bases used in Eq. (535) are 



M = 2> l e -l = 2' |eil = 7f' |e ~ l| = 7T (537) 
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The fact that the amplitude of the product is equal to the product of the amplitudes, as 
written in Eq. ( |514| ), can be demonstrated also by using a representation of the multiplication 
of the polar fourcomplex numbers by matrices, in which the polar fourcomplex number 
u = x + ay + j3z + jt is represented by the matrix 

x y z t 



(538) 



t x y z 
z t x y 
\y z t x ) 

The product u = x + ay + (3z + jt of the polar fourcomplex numbers u\ = x\ + ay\ + j3z\ + 
7*1, U2 = X2 + ay2 + (3z2 + 7*2, can be represented by the matrix multiplication 



A = At A 



1^2- 



It can be checked that the determinant det(A) of the matrix A is 



detA = v. 



(539) 



(540) 



The identity (514) is then a consequence of the fact the determinant of the product of matrices 
is equal to the product of the determinants of the factor matrices. 



5.3 The polar fourdimensional cosexponential functions 

The exponential function of a fourcomplex variable u and the addition theorem for the 
exponential function have been written in Eqs. ( |44"| ) and (|45|). If u = x + ay + (3z + 7*, then 
exp u can be calculated as exp u = exp x ■ exp(ay) • exp(/3z) • exp(7*). According to Eq. (|482|) , 



a 4m 


= l,a 4m+1 


= a, a 4m+2 = 0, a 4m+3 = 7, 




j3 2m 


= l,f3 2m+1 


= /?, 


(541) 






= 7,7 4m+2 = /3,7 4m+3 = a, 





where m is a natural number, so that exp(ay), exp(/3z) and exp(7z) can be written as 

exp(/3z) = coshz + /3sinhz, (542) 

and 

exp(ay) = 340(2/) + ag 4 i(y) + Pg&iv) + 7943 (y), (543) 
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exp(7t) = 340 (t) + 7341 (t) + f^m{t) + ag^(t), 



(544) 



where the four-dimensional cosexponential functions 340 , 341 , 342 , 343 are defined by the series 
340(2) = l + x 4 /4! + x 8 /8! + ---, (545) 



ff 4i (x) = x + x 5 /5! + x 9 /9\ + • • • , 
342(2) = x 2 /2\ + x 6 /6! + x 10 /10! + 



343 



(x) = x 3 /3! + x 7 /7! + x 11 /!!! + 



(546) 
(547) 
(548) 



The functions 340,342 are even, the functions 341,343 are odd, 



94o(-u) = 34o(«), 342(-«) = 342(«), 34i(-«) = -34i(«), 043(-«) = ~343(«)- (549) 



(550) 



It can be seen from Eqs. ( |545|) - (|548|) that 



340 + 341 + 342 + 343 = e , 340 - 341 + 342 - 343 = e" 



and 



340 - 342 = COSX, 341 - 343 = SIM, 



so that 



(340 + 341 + 342 + 343) (340 - 341 + 342 ~ 343) 



which can be also written as 



?40 - 342) 2 + (341 - 343) 2 



(551) 



1, (552) 



340-541 +342-f43- 2 (fi , 40f42-54l5'43)- 4 (5'4054l343 +34254l343 -34l5'40342-343540342) = 1.(553) 



The combination of terms in Eq. (553) is similar to that in Eq. ( |493[ ) . 

Addition theorems for the four-dimensional cosexponential functions can be obtained from 
the relation exp a(x + 3) = exp ax ■ exp ay, by substituting the expression of the exponentials 
as given in Eq. ( |543| ) , 



340 (x + y ) = 340 (2)340 (y) + 341 (2)343 (3) + 342 (2)342 (y) + m (2)341 (3) , 



341 (x + y) = 340 (2)341 (y) + 341 (2)340 (3) + 342 (2)343 (3) + 343 (2)342 (3) , 



(554) 



(555) 
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542 (x + y ) = 540 (SC)342 (y) + 541 (z)541 (y) + 542 (z)y40 (y ) + 543 (»543 (V ) , 



(556) 



543(2; + y) = 540(2?) 543 (y) + 541 (2^)542 (y) + 542(a;)54i(y) + 543(^)540(5)- 



For x = y the relations (|554] )- (|557D take the form 

540 (2ar) = 540(2;) + 542 0*0 + 2541(2)543 

541 (2x) = 2g 40 (x)g il (x) + 2542 (2)543 (2;), 

542 (2x) = 54 X (x) + 5I3 (x) + 254o(x)54 2 (x), 

543 (2x) = 254o(x)54 3 (x) + 254i(x)5 42 (x). 



For x = — y the relations ( |554|) -( |557| ) and ( |549| ) yield 

alo(x) + 542(2;) - 2541(^)543(2;) = l, 

541(2;) +543(2;) - 2540(2^)542(2;) = 0. 



From Eqs. ( 542[ )-( p44[ ) it can be shown that, for m integer, 
(cosh z + f3 sinh z) m = cosh mz + j3 sinh mz, 



(557) 

(558) 
(559) 
(560) 
(561) 

(562) 
(563) 

(564) 



and 



[54o(y)+«54i(y)+/3542(y)+7543(y)] m = 540 (™y)+«54i(ray)+/?542 (my) +7543(^5), (565) 

[540 (*) +7541 (t)+(3g42(t) + ag± 3 (t)] m = g m (mt)+ 7541 (mt)+ f3g4 2 (mt)+ ag 43 (mt). (566) 



Since 



(a + 7) 



2m r\2m— 1 



(1 + /3), (a + 7) 



2m+l i)2m 



2* n (a + 7 ), 



it can be shown from the definition of the exponential function, Eq. ( [44|) that 

1 + (3 a + 7 
exp(a + 7)x = ei H — cosh 2x H — sinh 2x. 



(567) 



(568) 



Substituting in the relation exp(a + j)x = expaxexp7x the expression of the exponentials 
from Eqs. ( 543| ), ( |544| ) and ( |568D yields 



2 2 2 2 1 + cosh 2x 

540 + 541 + 542 + 543 = o > 



(569) 



— 1 + cosh2x , r „. 
940942 + 941943 = t , (570) 

540941 + 940943 + 941942 + 942943 = £ sinh 2:r ' ( 571 ) 

where 940,941,942,943 are functions of x. 
Similarly, since 

(a - 7 ) 2m = (-l) m 2 2m " 1 (l - /?), (a - 7 ) 2m+1 = (-l) m 2 2m (a - 7), (572) 
it can be shown from the definition of the exponential function, Eq. ((44|) that 



1 + /3 

exp(a — 7)2; = — h e\ cos 2x + ei sin 2x. (573) 

Substituting in the relation exp(a— 7)2; = exp ax exp(-jx) the expression of the exponentials 
from Eqs. fljjjj ), and (07^) yields 

2 2,2 2 1 + COS2X (r.n A \ 
940 - 941 + 942 - 943 = 2 ' ^ > 

1 — cos2x ,.„, 

940942 - 941943 = ^ ' ( ' 

940941 - 940943 - 941942 + 942943 = ^ sin 2x, (576) 

where 940, 941, 942, 943 are functions of x. 

The expressions of the four-dimensional cosexponential functions are 

540 (x) = — (cosh x -\- cos x) , (577) 

541 (x) = - (sinh x + sin x), (578) 
942(2;) = -(coshx — cosx), (579) 
943(2;) = -(sinhx — sinx). (580) 

The graphs of these four-dimensional cosexponential functions are shown in Fig. 4. 

It can be checked that the cosexponential functions are solutions of the fourth-order 
differential equation 
d 4 C 

^ = C, (581) 

whose solutions are of the form Q{u) = Ag^u) + £941(1*) + Cg^.(u) + Dg^{u). It can also 
be checked that the derivatives of the cosexponential functions are related by 

^940 ^941 ^942 ^943 , cooN 

—j — = 943, —, — = 940, —j— = 941, — = 942- (582) 
aw aw aw aw 



5.4 The exponential and trigonometric forms of a polar four- 
complex number 

The polar fourcomplex numbers u = x + ay + f3z + 7* for which v + = x + y + z + t>0, v_ = 
x — y + z — t > can be written in the form 

x + ay + 0z + 7 t = e *i+<*2/i+/^i+7*i . ( 583 ) 

The expressions of xi,yi, z%,ti as functions of x,y,z,t can be obtained by developing 
e ayi ,e^ Zl and e 7 * 1 with the aid of Eqs. fl542|) -( |54"4]) , by multiplying these expressions and 
separating the hypercomplex components, and then substituting the expressions of the four- 
dimensional cosexponential functions, Eqs. (|577[) - (|580| ), 



x + y + z + t = e Xl+yi+Zl+tl , (584) 

x-z = e Xl - Zl cos(y 1 -t 1 ), (585) 

x-y + z-t = e *i-«i+*i-*i i (586) 

y-t = e x ^sm(y 1 -t 1 ). (587) 



It can be shown from Eqs. ( |584| )-( [587 ) that 



2 Z Z [1+ L 

where 

^ 2 = ( x + z ) 2 - (y + t) 2 = u+ > 0, v_ > 0. (589) 

The quantitites and w are determined by 

cos <f> = (x — z)/ fi + , sin <p = (y — t)//J,+ (590) 

and 

cosh a; = (x + z)/fi-, smhu = (y + i)//i_. (591) 

The explicit form of u> is 

1, x + y + z + t , 
u = -hx y . 592) 
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If u = u\u 2 , and = (xj + Zj) 2 — (yj + tj) 2 ,j = 1, 2, it can be checked with the aid of Eqs. 



(503) that 



P- = P1-P2- 



(593) 



Moreover, if cosh = (xj + Zj)//Jtj-, sinhujj = (yj + tj)/ fij-,j = 1, 2, it can be checked that 



ijj = U\ + U>2- 

The relation ( |594|) is a consequence of the identities 

(xtx 2 + z\z 2 + 2/1*2 + £12/2) + (^1^2 + ^1^2 + 2/12/2 + M2) 
= (xi + zi)(x 2 + 2 2 ) + (2/1 + £i)(2/2 + ^2), 



(594) 



(595) 



(x\y 2 + 2/1^2 + zit 2 + tiz 2 ) + (xit 2 + iix 2 + ^12/2 + yiz 2 ) 

= ( yi + tx)(x 2 + z 2 ) + (xi + Zl )(y 2 + t 2 ). (596) 

According to Eq. (|59C| ), is a cyclic variable, < <p < 2ir. As it has been assumed that 
v+ > 0, V- > 0, it follows that x + z > and x + z > \y + t\. The range of the variable u is 
—00 < uo < 00. The exponential form of the polar fourcomplex number u is then 



u = p exp 



-/3 In — + -a(w + ^) + ^7(^ 
2 /i+ 2 2 



where 



P = {P+P-) 



1/2 



(597) 



(598) 



If it = Ui«2, and pj = (fij+fij-) 1 / 2 , j = 1,2, then from Eqs. ( |509| ) and ( |593| ) it results that 



P = PlP2- 

It can be checked with the aid of Eq. fl533|) that 
1 + 13 



exp (ei( 



+ ei cos 4> + e± sin 0, 



(599) 



(600) 



which shows that e^ a ^l 2 is a periodic function of 4>, with period 27r. The modulus has the 
property that 



|uexp (ei4>)\ = \u\ 



(601) 



S3 



By introducing in Eq. ( |597| ) the polar angles 0+,8_ defined in Eqs. ( |515| ), the exponential 
form of the fourcomplex number u becomes 



u = p exp 



1 , m ^2 

- a + /3 + 7 In— — 
4 tan 8 , 



I( Q; - / 3 + 7 )l n _^- + e ~ 1 c 



(602) 



where < 9 + < vr/2,0 < 6L < 7r/2. The relation between the amplitude p, Eq. ( 598 ), and 
the distance d, Eq. ( |501| ), is, according to Eqs. (515) and ( |516| ), 



2 3 / 4 d 



(tan 6 + tan #_ 



1 1 
1 + - — — + 



tan 2 9+ tan 2 6. 



-1/2 



(603) 



Using the properties of the vectors e+,e_ written in Eq. (531), the first part of the expo- 

V2 



nential, Eq. (|602[) can be developed as 
exp 



1 72 1 

-(a + /3 + 7 )ln— - ---(a-/? + 7 )ln, 
4 tanfl+ 4 tan I 



^2 V2 
ex + e+- — h e_ 



1 \V4 
- tan 9 + tan 0_ 

2 / \ tant7_|_ tan l 
The fourcomplex number u, Eq. (|602| ) , can then be written as 



;/ = dV2 { 1 H + ' 



tan 2 6>+ tan 2 i 



-1/2 



ei + e H 



V2 



+ e_- 
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(604) 



exp (ei4>) , (605) 



tan 0+ ' tan #_ 
which is the trigonometric form of the fourcomplex number u. 

The polar angles 6 + ,9-, Eq. ( |515D , can be expressed in terms of the variables 9, A, Eq. 



(|523[) , as 

tan 9^ 
so that 



tan 9 cos A 



tan 9- 



1 



tan sin A 



1 1 

1 + o^T" + 



1 



tan 2 9j 



tan 2 cos 2 i 



(606) 



(607) 



The exponential form of the fourcomplex number u, written in terms of the amplitude p and 
of the angles 9,X,4> is 



u = p exp 



^(a + /3 + 7) ln(72tan6>cosA) - -{a- (3 + i) ln(72 tan (9 sin A) + en 



(608) 



where < A < ir/2. The trigonometric form of the fourcomplex number u, written in terms 
of the amplitude p and of the angles 9,X,(j) is 



u 



d\[2 [e\ cos 9 + e_|_ 72 sin 9 cos A + e_72 sin sin X^j exp (ei( 
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(609) 



If u = U1U2, it can be shown with the aid of the trigonometric form, Eq. (605), that the 
modulus of the product as a function of the polar angles is 



d 2 = U\d\ ( - + - „ ' 2n + 



1 1 1 



2 tan 2 9 1+ tan 2 <9 2 + tan 2 6»i_ tan 2 6 2 

1 + tan 2 9 1+ + tan 2 ) ( X + tan 2 2+ + tan 2 # 2 _) " (610) 
The modulus d of the product u\U2 can be expressed alternatively in terms of the angles 
9, A, (f) with the aid of the trigonometric form, Eq. ( |609| ) , as 

d 2 = Ad\d\ (— cos 2 6\ cos 2 6> 2 + sin 2 Q\ sin 2 62 cos 2 Ai cos 2 A 2 + sin 2 0\ sin 2 62 sin 2 Ai sin 2 A 2 

(611) 



5.5 Elementary functions of polar fourcomplex variables 

The logarithm ui of the polar fourcomplex number u, u\ = In it, can be defined as the 
solution of the equation 

u = e u \ (612) 

written explicitly previously in Eq. (|583| ) , for u\ as a function of u. From Eq. ( [597D it results 
that 

\nu = ln/9+ J^ln — + i«(o; + </>) + -7(0; - </>). (613) 

If the fourcomplex number u is written in terms of the amplitude p and of the angles 0+, 0_ , 0, 
the logarithm is 

1 \/2 1 \p2i 

In u = In p + - (a + (3 + 7) In — — - - (a - + 7) In — — - + ei0, (614) 
4 tanfl+ 4 tanfr_ 

where < 0+ < vr/2, < 9 + < ir/2. If the fourcomplex number u is written in terms of the 

amplitude p and of the angles 9, A, <j>, the logarithm is 

\nu = lnp+^(a + /3+7) ln(\/2 tan6> cos A) - i(a-/?+7) ln(\/2 tan 6> sin A) + ei0, (615) 

where 0<#<7r/2,0<A< tt/2. The logarithm is multivalued because of the term 
proportional to 0. It can be inferred from Eq. ( |614| ) that 

ln(iiiu 2 ) = lnui + ln«2, (616) 
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up to multiples of n(a — 7). If the expressions of p, and u in terms of x,y,z,t are 



introduced in Eq. (313), the logarithm of the polar fourcomplex number becomes 



In u = — ° ^ — — - ln(x+y+z+t)-\ — — ln(x— y+z— t)+e\ In /x + +ei^.(617) 

The power function u m can be defined for u + > 0, v _ > and real values of m as 

u m = e mlnu . (618) 
The power function is multivalued unless m is an integer. For integer m, it can be inferred 



from Eqs. ( 613 ) and ( |61Sj ) that 

( Ul u 2 ) m = u? v% . (619) 

Using the expression ( |617[ ) for Inu and the relations ( |531| )- (|534| ) it can be shown that 

u m = e + v+ +e-v™ + tfl (ei cos m(f) + e x sin m(f>) . (620) 

For integer m, the relation ( |620| ) is valid for any x, y, z, t. For natural m this relation can be 
written as 

u m = e+v™ + e.v™ + [ei(x - z) + h(y - t)\ m , (621) 

as can be shown with the aid of the relation 

ei cos mcf) + e\ sin mcj) = {e\ cos 4> + ei sin 0) m , (622) 

valid for natural m. For m = — 1 the relation ( |620| ) becomes 
1 _ 1 ^l + a + /3 + 7 l-a + /3-7 



x + ay + (3z + 7 i A\xA-y + z + t x — y + z — t 
l(l-/?)(x-z)-(a- 7 )(y-0 



(623) 



2 (x- z) 2 + (y- t) 2 

If m = 2, it can be checked that the right-hand side of Eq. (|621| ) is equal to (x + ay + /3z + 
7 t) 2 = x 2 + z 2 + 2yi + 2a(xy + zt) + /3(y 2 + t 2 + 2xz) + 2 7 (xt + yz). 

The trigonometric functions of the fourcomplex variable u and the addition theorems for 
these functions have been written in Eqs. (|90|)-(p^). The cosine and sine functions of the 
hypercomplex variables ay, f3z and jt can be expressed as 

cos ay = 540 - /%42, sin ay = ag Al - 7543, (624) 
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cos j3y = cosy, sin/??/ = /?siny, (625) 
cos jy = 540 - 0942, sin ^y = jg 41 - ag 4 3- (626) 
The cosine and sine functions of a polar fourcomplex number x + ay + (3z + 'yt can then be 



expressed in terms of elementary functions with the aid of the addition theorems Eqs. (92), 
(|93|) and of the expressions in Eqs. ( |624j) - (|626|) . 

The hyperbolic functions of the fourcomplex variable u and the addition theorems for 



these functions have been written in Eqs. (|97|)-( 100 ). The hyperbolic cosine and sine func- 



tions of the hypercomplex variables ay, (5z and jt can be expressed as 

cosh ay = g m + /3g 42 , sinh ay = aga + 7543, (627) 

cosh/3y = coshy, sinh/3y = /3sinhy, (628) 

cosh7y = 540 + (3g i2 , sinh -yy = 7g 41 + 0543. (629) 

The hyperbolic cosine and sine functions of a polar fourcomplex number x + ay + j3z + yt 
can then be expressed in terms of elementary functions with the aid of the addition theorems 
Eqs. (pG| ), ( |100D and of the expressions in Eqs. ( |627D -( |629|) . 



5.6 Power series of polar fourcomplex variables 

A polar fourcomplex series is an infinite sum of the form 

a + a\ + a 2 H h ai H , (630) 

where the coefficients ai are polar fourcomplex numbers. The convergence of the series ( |630| ) 
can be defined in terms of the convergence of its 4 real components. The convergence of 
a polar fourcomplex series can however be studied using polar fourcomplex variables. The 
main criterion for absolute convergence remains the comparison theorem, but this requires a 
number of inequalities which will be discussed further. 

The modulus of a polar fourcomplex number u = x + ay + /3z + jt can be defined as 

H = (x 2 + y 2 + 2 2 +* 2 ) 1/2 , (631) 
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so that according to Eq. (501) d = \u\. Since |x| < \u\, \y\ < \u\, \z\ < \u\, \t\ < \u\, a property 
of absolute convergence established via a comparison theorem based on the modulus of the 
series (630) will ensure the absolute convergence of each real component of that series. 

The modulus of the sum u\ + u 2 of the polar fourcomplex numbers ui,u 2 fulfils the 
inequality 



\U\\ - \u 2 \\ < \UI + U 2 \ < \ui\ + \U2\- 



(632) 



For the product the relation is 



l^i^l < 2\ui\ \u 2 \ 



(633) 



which replaces the relation of equality extant for regular complex numbers. The equality in 
Eq. fl633j ) takes place for x\ = y\ = z\ = ti, x 2 = J/2 = #2 = t 2 , or x\ = —y\ = z\ = —t\, x 2 = 
— 2/2 = z 2 = —t 2 . In particular 



\u 2 \ < 2(x 2 + y 2 + z 2 + t 2 ). 



(634) 



The inequality in Eq. fl633j) implies that 



lu'l < 2 l - 1 \u\ l . 



(635) 



From Eqs. ( 633 ) and (|635| ) it results that 



\au l \ < 2 l \a\\u\ l . 



(636) 



A power series of the polar fourcomplex variable u is a series of the form 



a + aiu + a 2 u 2 H h aiu 1 H . 



(637) 



Since 



1=0 



/=0 



a sufficient condition for the absolute convergence of this series is that 

Hm 2|a ; +i , IH < 1. 

Thus the series is absolutely convergent for 



(638) 



(639) 



\u\ < c , 



(640) 



where 



cq = lim 



\ a i\ 



(641) 



l— >oo 2 | 

The convergence of the series ( |637| ) can be also studied with the aid of the formula ( |621| 
which, for integer values of I, is valid for any x, y, z, t. If ai = ai x + aai y + (3ai z + 70/4, and 

= a ix + a iy + a lz + a it, 
M- = ai x — ai y + ai z - an, 
An = ai x — a iz, 
An = ai y — an, 



(642) 



it can be shown with the aid of relations ( |531| )-( |534| ) and ( 621 ) that the expression of the 



series (|637|) is 



00 

[Ai+v l + e + + Aj_v_e_ + (eiAn + (ei«i + e\V\) 1 



1=0 



(643) 



where the quantities v+ , V- have been defined in Eq. ( |495| ) , and the quantities v\ , «i have 



been defined in Eq. (|506|) . 

The sufficient conditions for the absolute convergence of the series in Eq. (|643j ) are that 



lim ^Y l \ V+l < 1, lim < 1. < 1. 

where the real and positive quantity Ai- > is given by 
Af = Af 1 + Al. 
Thus the series in Eq. (|643|) is absolutely convergent for 



\x + y + z + t\ < c+, \x — y + z — t\ < c_, /i+ < ci 



where 



Of = lim 



I A 



/+ 



\Ai 



1^00 \Ai + i 



lim 



-, ci = lim 



A, 



Z->oo Ai + i r 



(644) 



(645) 



(646) 



(647) 



The relations ( |646[) show that the region of convergence of the series ( 643 ) is a four-dimensional 
cylinder. It can be shown that cq = (1/2) min(c+, c_, ci), where min designates the smallest 
of the numbers in the argument of this function. Using the expression of \u\ in Eq. ( |508| ), it 
can be seen that the spherical region of convergence defined in Eqs. ( |640| ) , ( |641[) is included 
in the cylindrical region of convergence defined in Eqs. (|647| ). 
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5.7 Analytic functions of polar fourcomplex variables 

The fourcomplex function f(u) of the fourcomplex variable u has been expressed in Eq. ( |136| ) 
in terms of the real functions P(x, y, z, t), Q(x, y, z, t), R(x, y, z, t), S(x, y, z, t) of real variables 
x, y, z, t. The relations between the partial derivatives of the functions P, Q, R, S are obtained 
by setting succesively in Eq. (|l37l ) Ax — > 0, Ay = Az = At = 0; then Ay — > 0, Ax = Az = 
At = 0; then Az -> 0, Ax = Ay = At = 0; and finally At -> 0, Ax = Ay = Az = 0. The 
relations are 

FiC 

(648) 
(649) 
(650) 
(651) 

The relations ( |648[ )-( |65lD are analogous to the Riemann relations for the real and imag- 
inary components of a complex function. It can be shown from Eqs. d648|) -( |65l] ) that the 
component P is a solution of the equations 

Q2p Q2p ^ ^ Q2p Q2p ^ ^ 

dx 2 dz 2 ' dy 2 dt 2 

and the components Q, R, S are solutions of similar equations. As can be seen from Eqs. 

( 652| )-( |652[ ), the components P, Q, R, S of an analytic function of polar fourcomplex variable 

are harmonic with respect to the pairs of variables x,y and z,t. The component P is also a 

solution of the mixed-derivative equations 

c^P_a^P d 2 P _ d 2 P #P_9^P d 2 P _ d 2 P 
dx 2 dydt ' dy 2 dxdz ' dz 2 dydt ' dt 2 dxdz ' 

and the components Q, R, S are solutions of similar equations. The component P is also a 

solution of the mixed-derivative equations 

Q2p Q2p Q2p Q2p 



dP 


_ dQ _ 


dR 


dS 


dx 


dy 


dz 


~ 'dt 1 


dQ 


dR 


dS 


dP 


dx 


dy 


dz 




OR 


dS 


dP 


dQ 


dx 


dy 


dz 


~ dt ' 


OS 


dP 


dQ 


dR 


dx 


dy 


dz 


~ ~dt' 



dxdy dzdt' dxdt dydz" 1 
and the components Q,R, S are solutions of similar equations. 



(654) 
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5.8 Integrals of functions of polar fourcomplex variables 

The singularities of polar fourcomplex functions arise from terms of the form l/(u — ito)"\ 
with m > 0. Functions containing such terms are singular not only at u = Uq, but also at all 
points of the two-dimensional hyperplanes passing through no an d which are parallel to the 
nodal hyperplanes. 

The integral of a polar fourcomplex function between two points A, B along a path 
situated in a region free of singularities is independent of path, which means that the integral 
of an analytic function along a loop situated in a region free from singularities is zero, 

f(u)du = 0, (655) 

r 

where it is supposed that a surface S spanning the closed loop T is not intersected by any of 
the hyperplanes associated with the singularities of the function f(u). Using the expression, 
Eq. (|136j ) for f(u) and the fact that du = dx + ady + (3dz + jdt, the explicit form of the 
integral in Eq. ( |655| ) is 

j> f(u)du = j> [(Pdx + Sdy + Rdz + Qdt) + a(Qdx + Pdy + Sdz + Rdt) 

+(3{Rdx + Qdy + Pdz + Sdt) + j(Sdx + Rdy + Qdz + Pdt)\ . (656) 

If the functions P,Q,R,S are regular on a surface S spanning the loop T, the integral along 
the loop r can be transformed with the aid of the theorem of Stokes in an integral over the 
surface S of terms of the form dP/dy — dS/dx, dP/dz — dR/dx, dP/dt — dQ/dx, dR/dy — 
dS/dz, dS/dt — dQ/dy, dR/dt — dQ/dz and of similar terms arising from the a,f3 and 7 
components, which are equal to zero by Eqs. ( 648| )-( |65TD , and this proves Eq. fl655| ). 



The integral of the function (u — uo) m on a closed loop T is equal to zero for m a positive 
or negative integer not equal to -1, 

<j) (u — uo) m du = 0, m integer, m/-l, (657) 

This is due to the fact that J(u — uo) m du = (u — uo) m+1 j(ra + 1), and to the fact that the 
function (u — uo) m+1 is singlevalued for m an integer. 

The integral fp du/(u — uq) can be calculated using the exponential form (|602| ). 



u — uq = p exp 



ha + (3 + 7) In — ^- - ha - + 7) In — ^- + en 
4 tanfl+ 4 tanfl_ 



(658) 
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so that 



''" <,P + i(a + /3 + 7)dln-^-i(a- / a + 7 )dln-^ + ei#- (659) 



u — uq p 4 tan 0+ 4 tan 0. 

Since p, tan # + and tan are singlevalued variables, it follows that fpdp/p = 0,^dln y/2/ tan i 
0, and fpdlny/2/ tan# + = 0. On the other hand, is a cyclic variables, so that it may give 
a contribution to the integral around the closed loop T. The result of the integrations will 
be given in the rotated system of coordinates 

£ = ^=0 - z), v = -j=(y - t), t = ~(x + y + z + t), v = ^(x - y + z - t). (660) 

Thus, if Cm is a circle of radius r parallel to the £Ov plane, and the projection of the center 
of this circle on the £Ov plane coincides with the projection of the point uq on this plane, 
the points of the circle Cii are described according to Eqs. (|506| )- (l507j ) by the equations 

£ = Co + r cos <f>, v = v + r sin r = r , ( = Co, (661) 

where uq = xo + ayo + (5zq + 7*0, and £0, f o> r o> Co ar e calculated from xo, yo> -^o, *o according 
to Eqs. fl66CD . 
Then 

du 

2vrei. (662) 



/C|| U — Uq 

The expression of §pdu/(u — uq) can be written with the aid of the functional int(M, C) 
defined in Eq. ( |153[) as 
f du 

f = 27rei int(w £„,r£„), (663) 

Jr u — u 

where uq^ v and T^ v are respectively the projections of the point uq and of the loop T on the 
plane £v. 

If f(u) is an analytic polar fourcomplex function which can be expanded in a series as 



written in Eq. (131), and the expansion holds on the curve T and on a surface spanning T, 
then from Eqs. ( |657| ) and ( |663| ) it follows that 
f f(u)du 



lr u-u 



27rei mt(u oiv ,T iv )f(u ), (664) 
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where T^ v is the projection of the curve T on the plane £v, as shown in Fig. 5. Substituting 



in the right-hand side of Eq. (664) the expression of f(u) in terms of the real components 
P,Q,R,S, Eq. (pe|), yields 

r fMdu = 7r[{(3 _ 1)(Q _ 5) + (q _ 7)(P _ int (n r ) (665) 

Jr u — u 

where P, Q, -R, S are the values of the components of / at u = uq. If the integral is written 



as 



= i + a i a + /3/ + 7 J (666) 

IT U - Uq 

it results from Eq. ( |665| ) that 

I + I a + Ip + 1 7 = 0. (667) 
If f(u) can be expanded as written in Eq. ( |131| ) on T and on a surface spanning T, then 



from Eqs. (|657|) and (|663| ) it also results that 

f(u)du 7T . „ . „/_'i , 

, \ m+1 = 2—~ei int ( UQiv ,T iv ) n , 668 
r (u — no) + m! 

where the fact has been used that the derivative f^ m '(uo) of order m of /(u) at u = uq is 
related to the expansion coefficient in Eq. (|131| ) according to Eq. ( |135| ). 

If a function /(n) is expanded in positive and negative powers of u — Uj, where Uj are polar 
fourcomplex constants, j being an index, the integral of / on a closed loop T is determined 
by the terms in the expansion of / which are of the form <Xj/{u — Uj), 



/(„) = . + + (669) 

~ U — Uj 

3 J 

Then the integral of / on a closed loop T is 

f f(u)du = 27rei int(uj£ v , r^ v )aj. (670) 

5.9 Factorization of polar fourcomplex polynomials 

A polynomial of degree m of the polar fourcomplex variable u = x + ay + j3z + jt has the 
form 

P m (u) = u m + a lU m ~ l + ■■■ + a m -iu + am, (671) 
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where the constants are in general polar fourcomplex numbers. 



If a m = a mx + aa my + (3a mz + 7a m t, and with the notations of Eqs. (495) and (p4z 
applied for 0, 1, ■ ■ ■ , m , the polynomial P m (u) can be written as 



v rn + AlV m-l + . . . + Am _ lV+ + At 



e+ + 



v rn + Al , v rn-l + ... + A ^_ iV _ + A ^ 



+ 



IIL 

(eivi + eivi) m + ( e lAl + ) ( e i' u i + 
" i 



,m—l 



(672) 



where the constants ^4;+, ^4;_ , A^, An are real numbers. The polynomial of degree m in 
{e\V\ + &\V\) can always be written as a product of linear factors of the form [ei(ui — vu) + 
— where the constants vu,vn are real. The two polynomials of degree m with 
real coefficients in Eq. ( |672| ) which are multiplied by e+ and e_ can be written as a product 
of linear or quadratic factors with real coefficients, or as a product of linear factors which, if 
imaginary, appear always in complex conjugate pairs. Using the latter form for the simplicity 
of notations, the polynomial P m can be written as 



p m = ~[[(v+ - s/+)e+ + Y[(v - s;-)e_ + [ei(«i - v u ) + e x {vi - vu)} 



(673) 



l=i i=i i=i 

where the quantities si + appear always in complex conjugate pairs, and the same is true for 

the quantities s/_. Due to the properties in Eqs. fl53lD -( |534| ), the polynomial P m (u) can be 

written as a product of factors of the form 



P m (u) = [(«+ - s /+ )e+ + (t 
i=i 



■ s/_) e_ + (ei(«i - vu) + ei(vi - vu))] . 
These relations can be written with the aid of Eq. ( ^35[) as 



P m (u) = Y[(u - u p ), 
P =i 



where 



Up = s p+ e + + Sp-e- + e\V\ v + e\V\ v ,-p = 1, m. 



(674) 



(675) 



(676) 



The roots s p+ , s p _, vi p e% + vi p e\ of the corresponding polynomials in Eq. (373) may be 
ordered arbitrarily. This means that Eq. ( |676| ) gives sets of m roots u\, ...,u m of the poly- 
nomial P m (u), corresponding to the various ways in which the roots s p+ , s p _, v\ p e\ + v Xp e\ 
are ordered according to p in each group. Thus, while the hypercomplex components in 
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Eq. (672) taken separately have unique factorizations, the polynomial P m (u) can be written 
in many different ways as a product of linear factors. The result of the polar fourcomplex 
integration, Eq. fl670|), is however unique. 



If P(u) 



1, the factorization in Eq. (|675|) is u 2 — 1 



u\){u — u 2 ), where 



+e + ± e_ ± ei,U2 = — iti, so that there are 4 distinct factorizations of u 2 — 1, 
u 2 - 1 = (u - l)(u + 1), 



l = (n-/3)(u + /3), 



u 2 - 1 



2 

-1+Q+/3+7 
2 



(677) 



It can be checked that {±e+ ± e_ ± e\Y 



e+ + e_ + ei 



5.10 Representation of polar fourcomplex numbers by irre- 
ducible matrices 



If T is the unitary matrix, 



/ i 

2 

1 
2 



^2 U v 7 ! 



V 



^2 



(678) 



TUT- 



\ 



(679) 



it can be shown that the matrix TUT 1 has the form 
( x+y+z+t ^ 

x-y+z-t 

J 

where U is the matrix in Eq. ( |533| ) used to represent the polar fourcomplex number u. In 
Eq. (679), V\ is the matrix 
x — z y — t 
-y + t x — z 

The relations between the variables x+y+z+t, x—y+z—t,x — z,y—t for the multiplication of 



Vt 



(680) 



polar fourcomplex numbers have been written in Eqs. ( 504 ),( p0q ), ( ^12| ), ( |513| ). The matrix 
TUT -1 provides an irreducible representation ^ of the polar fourcomplex number u in terms 
of matrices with real coefficients. 
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6 Conclusions 



In the case of the circular fourcomplex numbers, the operations of addition and multiplication 
have a simple geometric interpretation based on the amplitudes p, the azimuthal angles (ft, x, 
and the planar angles tft. Exponential form exist for the circular fourcomplex numbers, 
involving the variables v, (ft, \ and tft- The circular fourcomplex functions defined by series of 
powers are analytic, and the partial derivatives of the components of the circular fourcomplex 
functions are closely related. The integrals of circular fourcomplex functions are independent 
of path in regions where the functions are regular. The fact that the exponential forms of 
the circular fourcomplex numbers depend on the cyclic variables (ft, % leads to the concept 
of pole and residue for integrals on closed paths. The polynomials of circular fourcomplex 
variables can always be written as products of linear factors. 

In the case of the hyperbolic fourcomplex numbers, an exponential form exists for the 
hyperbolic fourcomplex numbers, involving the amplitude p and the arguments y\,z\,t\. The 
hyperbolic fourcomplex functions defined by series of powers are analytic, and the partial 
derivatives of the components of the hyperbolic fourcomplex functions are closely related. 
The integrals of hyperbolic fourcomplex functions are independent of path in regions where 
the functions are regular. The polynomials of tricomplex variables can be written as products 
of linear or quadratic factors. 

In the case of the planar fourcomplex numbers, the operations of addition and multipli- 
cation of the planar fourcomplex numbers introduced in this work have a geometric inter- 
pretation based on the amplitude p, the azimuthal angles (ft, x, and the planar angle tft. An 
exponential form exists for the planar fourcomplex numbers, involving the variables p, (ft, x 
and ift. The planar fourcomplex functions defined by series of powers are analytic, and the 
partial derivatives of the components of the planar fourcomplex functions are closely related. 
The integrals of planar fourcomplex functions are independent of path in regions where the 
functions are regular. The fact that the exponential form of the planar fourcomplex numbers 
depends on the cyclic variables (ft, x leads to the concept of pole and residue for integrals 
on closed paths. The polynomials of planar fourcomplex variables can always be written as 
products of linear or quadratic factors. 

In the case of the polar fourcomplex numbers, the operations of addition and multipli- 
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cation have a geometric interpretation based on the amplitude a, the polar angles 9+,9- 
and the azimuthal angle 0. An exponential form exists for the polar fourcomplex numbers, 
involving the variables a, (ft and the angles The polar fourcomplex functions defined 

by series of powers are analytic, and the partial derivatives of the components of the polar 
fourcomplex functions are closely related. The integrals of polar fourcomplex functions are 
independent of path in regions where the functions are regular. The fact that the exponential 
form of the polar fourcomplex numbers depends on the cyclic variable <j> leads to the concept 
of pole and residue for integrals on closed paths. The polynomials of polar fourcomplex 
variables can be written as products of linear or quadratic factors. 
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FIGURE CAPTIONS 



Fig. 1. Azimuthal angles </>, \ an d planar angle ip of the fourcomplex number x + ay + 
j3z + 7t, represented by the point A, situated at a distance d from the origin O. 

Fig. 2. Integration path V and the pole uq, and their projections r^„,r T ^ and uq^ v ,uq ti ^ 
on the planes £v and respectively t£. 

Fig. 3. The planar fourdimensional cosexponential functions /40, /41, /42, /43- 

Fig. 4. The polar fourdimensional cosexponential functions g^o, 541, 542, 543. 

Fig. 5. Integration path T and the pole no, and their projections T^ v and u ^ v on the 
plane 
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Fig. 5 



